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ABSTRACT 


A  theory  of  photodesorption  via  resonant  laser  excitation  of 
internal  adsorbate  molecular  vibrations  is  presented.  It  is  applied 
to  physisorbed  molecules  at  low  surface  coverage  on  both  dielectric 
and  perfect  conductor  substrates. 

Effects  of  the  laser,  phonons,  and  tunneling  are  included  in  a 
calculation  of  microscopic  transition  rates  among  the  states  of  the 
physisorp tion  surface  bond.  Derivation  of  the  macroscopic  desorption 
rate  is  achieved  by  numerical  solution  of  a  master  rate  equation  incor¬ 
porating  the  microscopic  phenomena.  Desorption  is  found  to  be  both 
temperature  and  laser  intensity  dependent,  from  low  temperatures  up  to 
the  thermal  desorption  regime,  and  for  intensities  up  to  saturation. 

The  influence  of  the  system  parameters  -  including  the  molecular  vibra¬ 
tional  frequency  and  the  depth  and  range  of  the  physisorp tion  potential- 
on  desorption  is  also  elucidated.  Finally,  results  presented  for  the 
CH^F/NaCl  and  GO/Cu  systems  are  compared  with  experiments  in  the 
li terature . 
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§1.  INTRODUCTION 


The  kinetics  of  gas-solid  surface  interactions  has  been  a  subject 
of  escalating  interest  -  both  experimental  and  theoretical  -  in  recent 
years.  Chemical  reactions  catalyzed  at  solid  surfaces  are  of  funda¬ 
mental  importance  -  to  industrial  applications,  for  example  -  and  their 
kinetic  theory  is  under  current  development.  The  processes  of  adsorp¬ 
tion  and  desorption  of  gas  molecules  at  a  solid  surface  are  an  important 
part  of  complex  catalytic  reactions. 

Adsorption  is  the  process  whereby  gas  phase  particles  enter  into 
an  energetically  favourable  configuration  at  a  solid  surface.  (Desorp¬ 
tion  occurs  when  gas  particles  leave  this  configuration.)  If  the  bond 
formed  between  the  gas  and  surface  molecules  is  strong,  with  binding 
energy  typically  greater  than  1  eV,  then  the  reaction  is  called  chemi¬ 
sorption.  In  this  case  extensive  "electron  sharing"  occurs  between  gas 
and  solid  atoms.  While  this  reaction  is  of  prime  importance  to  cataly¬ 
sis,  we  consider  instead  the  weaker  gas-surface  interaction  of  physi- 
sorption.  A  gas  molecule  is  said  to  be  physisorbed  at  a  surface  if  the 
bond  is  weak  -  typically  with  binding  energy  <1  eV.  In  this  case 
electron  sharing  is  minimal,  and  the  gas  molecules  are  chemically 
unaltered  after  undergoing  an  adsorption  and  desorption  process.  We 
will  focus  our  attention  on  the  desorption  of  physisorbed  polyatomic 
molecules . 

Desorption  may  proceed  via  many  different  channels;  two  of  major 
current  interest  are  summarized  below: 
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i)  The  thermal  motion  of  the  solid  surface  may  couple  to  the  gas- 

solid  bond.  See  Fig.  la.  By  absorption  of  phonons,  the  surface  bond 

may  be  excited  to  a  continuum  state,  in  which  case  desorption  occurs. 

This  process,  called  thermal  desorption,  has  been  analyzed,  using  a 

1  2 

full  quantum  statistical  theory,  by  Gortel  and  Kreuzer  *  .  Two  common 

methods  for  initiating  the  reaction  are  1)  reducing  the  gas  pressure 

at  the  surface,  and  2)  heating  the  solid  either  slowly  ("programmed 

thermal  desorption")  or  rapidly  ("flash  desorption") . 

Experiments  show  that  the  single  gas  particle  desorption  rate  may, 

3 

in  a  certain  regime,  be  parameterized  by  an  Arrhenius- type  formula  : 

-Ed/kT 


-1 

t ,  =  v  e 

a 


(1) 


where  t^  is  the  characteristic  desorption  time,  T  is  the  solid  temper¬ 
ature,  is  the  heat  of  adsorption,  and  v  is  a  prefactor  that  is 

classically  approximated  by  the  vibrational  frequency  of  the  particle 
in  the  surface  bond.  The  more  general  theory  of  Gortel  and  Kreuzer^ 
predicts  desorption  via  a  cascade  of  step  by  step  phonon-assisted 
transitions  up  through  the  surface  bond  states  into  the  continuum. 

Their  results,  which  agree  well  with  experiments,  also  delineate  the 
regime  in  which  the  simple  formula  given  above  may  be  applied, 
ii)  Another  means  of  stimulating  desorption  is  to  irradiate  the  sur¬ 
face  with  high  intensity  laser  light  -  hence  photodesorption.  The 
radiation  energy  may  couple  to  the  adsorbate-solid  system  either  1)  by 

directly  heating  the  solid  via  non- resonant  absorption  of  radiation 

4 

from  a  large  laser  pulse  ,  2)  by  direct  excitation  of  the  adsorbate- 
5  6 

surface  bond  ’  ,  or  3)  by  resonant  excitation  of  an  internal  vibrational 
mode  of  an  adsorbate  molecule^  (See  Figs.  lb,c,d,  respectively.) 
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Fig.  1.  Possible  channels  to  desorption. 
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We  will  only  consider  systems  at  low  surface  coverage,  and  have  there¬ 
fore  ignored  the  possibility  of  interactions  between  the  adsorbate 
molecules  themselves. 

Reaction  mechanism  (1)  above  proceeds  via  thermal  desorption  from 
a  suddenly  heated  surface,  and  has  been  termed  "laser  flash  desorption". 
On  the  other  hand,  the  surface  bond  itself  absorbs  the  laser  radiation 
in  channel  (2).  Finally,  reaction  channel  (3)  permits  the  resonant 
absorption  of  light  by  the  adsorbate  molecule  itself,  with  subsequent 
tunneling  processes  leading  to  desorption.  It  is  anticipated  that 
mechanism  (3)  will  generally  be  more  efficient  than  (2),  because  anhar- 
monicity  of  the  surface  potential  makes  light  absorption  in  (2)  largely 
non-resonant.  This  is  explained  in  greater  detail  in  Section  3.3. 

We  concentrate  on  channel  (3)  above.  Since  different  species  of 
adsorbate  molecules  have  different  radiation  absorption  spectra, 
mechanism  (3)  may  possibly  be  used  to  selectively  excite  the  particu¬ 
lar  species  we  wish  to  desorb,  by  tuning  the  laser  frequency  appro¬ 
priately.  The  possibility  then  exists  of  using  this  mechanism  to 
purify  and  separate  substances,  or  to  determine  the  individual  mole¬ 
cular  populations  of  a  multicomponent  adsorbate. 

Some  basic  details  of  the  gas-solid  system  that  we  consider  are 
now  presented.  It  is  reminded  that  we  deal  only  with  a  physisorbed, 
low  coverage  adsorbate  at  very  low  gas  pressures.  We  assume  further 
that  the  adsorbate  molecules  reside  at  localized  sites  and  are  not 
free  to  move  laterally  on  the  surface.  The  solid  is  treated  as  a  con¬ 
tinuum,  which  interacts  with  the  adsorbate  via  its  surface  only.  Our 
adsorbate  consists  of  diatomic  molecules  adsorbed  end  on  -  that  is, 
with  one  atom  much  closer  to  the  surface  than  the  other.  We  may  model 
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some  polyatomic  molecules  too,  by  replacing  the  outer  "atom"  above  with 
a  complex  (e.g.  F-CH^) .  The  rotational  motion  of  gas  molecules,  whether 
adsorbed  or  free,  is  ignored.  By  ignoring  surface  heating  effects  via 
channel  (1)  above,  the  choice  of  solid  is  restricted  to  either  a  di¬ 
electric  which  is  transparent  at  laser  frequencies,  or  to  a  very  good 

conductor.  (In  the  latter  case  we  must  also  be  careful  to  not  use  too 

8b 

high  laser  intensities  .)  Finally,  our  model  process  for  resonant 
photodesorption  requires  a  laser  that  can  be  tuned  to  the  infrared 
spectral  line  of  the  appropriate  vibrational  mode  of  an  adsorbate  mole¬ 
cule  . 

Within  the  bounds  of  this  model,  we  will  calculate,  in  a  quantum 
mechanical  formalism,  microscopic  transition  rates  from  which  the 
macroscopic  desorption  phenomena  can  be  derived;  see  Section  2.  The 
phenomenological  characteristics  of  our  model  are  studied  in  Section  3, 
where  a  comparison  is  made  with  presently  available  experimental 
results  and  with  other  theories  in  the  literature.  Concluding  remarks 


are  found  in  Section  4. 


•: 


§2.  CALCULATION  OF  THE  PHOTODESORPTION  PATES 


§2.1  Summary  of  Program 

We  want  to  treat  the  interaction  process  of  Fig.  2  quantum  mechan¬ 
ically,  and  develop  rate  equations  to  describe  the  desorption  process. 


e.m.  field  * — *  internal  molecular 
(photons)  vibrations 


tunne ling 


Lsolid  vib rations-! 
(phonons) 


surface  bond 
vib ra tions 


Fig.  2.  The  interaction  process. 


Suppose  the  surface  bond  is  described  by  a  potential  well  (see 
Fig.  3a)  which  has  bound  (adsorbed)  states  labelled  by  the  quantum 
number  j,  and  a  continuum  of  free  states,  each  denoted  q.  (We  will 
use  a  Morse  potential  to  model  the  surface  bond.)  Further,  assume  that 
internal  vibrations  of  the  adsorbed  gas  molecule  are  approximately  sim¬ 
ple  harmonic,  giving  an  energy  spectrum  of  equally  spaced  levels  denoted 
by  the  vibrational  state  index  v.  The  double-index  (j,v)  fully  des¬ 
cribes  the  state  of  an  adsorbate  molecule  in  our  theory.  The  inter¬ 
action  process  in  Fig.  2  causes  transitions  between  the  states  (j,v). 

If  a  gas  molecule,  previously  adsorbed  in  the  state  (jf,vf),  is  trans¬ 
ferred  to  a  free  state  (q,v),  then  the  surface  bond  is  broken  and  the 
molecule  leaves  the  surface  with  some  translational  velocity.  This 
process  is  called  desorption.  Note  that  an  adsorbate  molecule  in  state 
(j,v)  may  have  total  energy  greater  than  zero,  if  the  positive  contri¬ 
bution  due  to  the  internal  vibrational  state  (v)  exceeds  the  negative 
portion  due  to  the  bound  surface  state  (j).  Whether  the  molecule  is 
adsorbed  or  free  depends  on  the  surface  bond  index  (j  or  q)  only.  As 
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7 


Fig.  3a.  Energy  diagram  of  gas  molecule  states. 


Fig.  3b.  Two  possible  paths  to  desorption 
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an  example,  two  possible  transition  paths  to  desorption  are  shown  gra¬ 
phically  in  Fig.  3b.  Notice  that  the  path  to  desorption  may  be  very 
complex,  involving  many  transitions  along  the  way. 

We  assume  that  the  dynamics  of  the  adsorbate  molecules  is  a  sta¬ 
tionary  Markov  process.  Defining  the  state  occupation  numbers  nY(t) 

•3 

we  describe  the  time-evolution  of  our  system  by  the  master  equation 


vv’  v* 


J  (v'.J'Wv.j)  33  J 


R  I  (t)  - 


i  rv!v+  i  rv:v 

3  3  vVv  C3 


nT(t)  . 


(2) 


where  rTT,  is  the  single  particle  transition  rate  from  (j',v')  to  (j,v), 
w* 

R  is  the  single  particle  transition  rate  from  (j',vf)  to  all  con- 

tinuum  states  (q)  with  vibrational  index  v.  (rTT,  and  R^T,  are  given 

J  J  c3 

by  our  microscopic  quantum  mechanical  theory  described  later  in  this 
section. ) 


We  assume  that,  initially,  the  adsorbate  is  already  present  on  the 
surface  and  that  no  further  adsorption  occurs  during  the  desorption  pro¬ 
cess.  That  is,  there  are  no  transitions  (q,  v* )  -+  ( j ,v) .  This  is  reason¬ 
able  since  photodesorption  experiments  are  normally  performed  in  ultra 
high  vacuum  conditions;  there  are  exceedingly  few  free  gas  molecules 
available  for  adsorption. 

In  postulating  the  rate  equations  (2),  we  assume  that  the  thermal 
equilibrium  in  the  solid  and  the  steady  state  conditions  of  the  laser 
are  maintained  throughout  the  desorption  process.  The  role  of  equations 
(2)  in  a  full  quantum  statistical  formalism  is  elucidated  in  reference  1. 

The  desorption  rate  is  now  just  the  rate  of  depopulation  of  the 


bound  states  (j),  and  hence: 


_d_ 

dt 


(3) 


Des.rate  =  N  t 

d 


where 


I  »T(t)-  I  i!(t) 

(j,v)  J  (j,v)  J 

V1  v 

t  is  the  characteristic  desorption  time  and  N=  /,  n.(0). 

(j»v)  J 

In  principle,  the  desorption  rate  may  be  determined  by  solving 
the  coupled  system  of  rate  equations  (2)  for  specified  initial  condi¬ 
tions.  We  rewrite  (2)  in  more  compact  form  as 


->■  -> 
n(t)  =  Rn(t) 


(4) 


where  n(t)  is  a  column  vector  of  occupation  numbers  n.(t)  and  R  is  the 
rate  matrix,  with  element 


R .  .  t  i  — 


rTTi 

33 


■  I  rYL'Y  -Jr. 

(v",j")  J  J  V'  CJ 


v"v 


for  (j,v)  ^(j'.v1) 
for  diagonal  elements 


(We  include  only  bound  states  (j)  in  the  rate  equations  since  we  ignore 
the  possibility  of  transitions  returning  from  the  continuum.) 

Solving  (4)  involves  a  complex  numerical  calculation,  but  we  may  esti- 

<->■ 

mate  the  desorption  time  by  studying  the  eigenvalues  of  matrix  R,  as 
described  in  Section  2.5.^ 

w1  w 1 

It  remains  to  calculate  the  transition  rates  R...  and  R  using 

33  cj’ 

the  interaction  channel  described  in  Fig.  2.  To  include  the  laser 
radiation  correctly,  we  must  first  determine  the  electromagnetic  field 
near  the  solid  surface  in  a  suitable  quantized  form;  see  Section  2.2. 
The  Hamiltonian  of  the  gas-solid-laser  system,  including  appropriate 
interaction  terms,  is  constructed  in  Section  2.3.  A  straightforward 
application  of  time-dependent  perturbation  theory  then  gives  the  tran¬ 


sition  rates  in  Section  2.4. 


§2.2  Quantization  of  the  Laser  Radiation 


a)  Introduc tion 

We  shall  calculate  the  quantized  electromagnetic  field  for  the 
special  case  in  which  the  solid  is  a  lossless  non-magnetic  dielectric, 
and  then  at  the  end  show  how  one  can  solve  the  problem  if  the  solid 

•k 

were  a  perfect  conductor  instead. 

Suppose  the  dielectric  filling  the  half  space  x<  0  is  described 
by  a  uniform,  scalar  permitivity  e,  and  the  x>0  region  is  the  vacuum. 
(See  Fig.  4a) .  We  ignore  the  effect  of  gas  molecules  on  the  electro¬ 
magnetic  field,  since  the  molecular  density  has  already  been  assumed 
to  be  small,  both  on  the  solid  surface,  and  in  the  free  gaseous  state. 


Fig.  4a.  The  semi-infinite  solid. 


We  can,  in  fact,  obtain  the  perfect  conductor  results  from  those 
for  the  dielectric  case  by  letting  the  refractive  index  (n  )  become 
infinite . 


■ 
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To  determine  the  electromagnetic  field  near  the  solid  surface,  one 
solves  Maxwell’s  equations  for  the  appropriate  geometry  and  boundary 


conditions.  These  equations  are: 

V-E  =0,  V-S=0, 

with 

f  xn,  n2(x)  Is  •  n,  B  continuous  at  x  =  0 


-±  9  it 

v  xE=~’ 


„  2 .  .  3E 

V  x  B  =  n  (x)  e  \x  — 
o  o  dt 


(5a) 


(5b) 


where 


n(x)  = 


1  ;  x  >  0 


[nr  ;  x  <0 


FT 

;  n  =  I— 


r  v 


is  the  refractive  index  of  the 
dielectric. 


The  solutions  to  (5)  must  be  consistent  with  incident  plane  waves  from 
the  x>0  side.  Combining  (5a)  gives  the  wave  equations 


V2u  -  — ^ — ^-4  =  0 


v2(x)  9t2 


(6) 


Q 

where  u  is  a  component  of  E  or  B,  and  v(x)  =  — —  is  identified  as  the 


n(x) 


wave  velocity. 


-+  ,-y 


Introducing  the  vector  potential,  A(x,t),  in  the  Coulomb  gauge  in 
the  usual  manner,  we  have 


|  _  JA 

E  __9t  * 


B  =  Vx A  , 


V-A  =  0 


(7) 


Substituting  (7)  into  (5a, b) ,  the  problem  is  rewritten  as 
1 


02-> 

V  A  - 


2-v 

9  A 


2  2 

v  (x)  9t 


=  0 


(8a) 


where 


~y  ~y  -y  2 ,  s  ~y  „  n 

n  xa  ,  V  xA  ,  n  (x)  A*n  are  continuous  at  x  =  0  . 


(8b) 


For  the  quantization  procedure  (described  formally  in  (c))  we  will 
need  to  expand  the  electromagnetic  potential  A(x,t)  in  terms  of  a  set  of 
orthogonal  modes  {U-»-  (x)  }  which  independently  satisfy  equations  (8)  : 

he 


« 


- 
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A(x,  t)  =  T  $  U  (x)  e  C 

it,3  1  &  ] 


leg 


+  hermitian  conjugate 


(9) 


where  k  and  3  are  the  incident  wave  vector  and  polarization  of  the 

mode,  is  a  normalization  constant,  C->  is  the  mode  amplitude,  that 

k3  k3 

will  become  the  annihilation  operator  after  quantization  and  £1*  is  the 


k3 


oscillation  frequency. 

In  Section  (b)  we  construct  a  suitable  set  of  modes  {U->  (x)  }  and 


k3 


present  the  appropriate  orthogonality  relations  needed  in  the  quanti¬ 
zation  formalism  described  in  Section  (c). 

b)  Orthogonal  modes 

To  find  solutions  of  equations  (8)  appropriate  to  the  laser  radia¬ 
tion,  we  consider  a  plane  wave  incident  on  the  dielectric  at  angle  0, 
as  shown  in  Fig.  4b.  Instead  of  solving  (8)  directly,  we  determine  the 


Fig.  4b.  Ray  diagram  of  laser  radiation. 


' The  sum  in  equation  (9)  is 
infinite  size,  the  mode  wave 


in  (9)  is  symbolic  for  L  die.  Since  the  system  is  of 
mode  wave  number  (k)  i^  continuous. 
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reflected  and  transmitted  waves  using  Snell's  law  and  Fresnel's  equa¬ 
tions11.  The  wave  vectors  of  the  incident,  reflected  and  transmitted 
waves  are  respectively: 


k.  =  KK  +  £n  ,  k 
1  i 


KK  -  Cn  ,  it  =  KK  -  bn 


(10) 


a  y  *  ^  - 

where  £  <0;  K,n  are  shown  in  Fig.  4b,  and  b  =n  k.  cos  0  = 

°  r 1  1 1  t 

/(nr-l)  K  +  n^  £  by  Snell's  law.  We  use  the  coordinate  system  x  =  R+ nx 

where  ll*n  =0.  Let  the  incident  radiation  be  described  by  (^S  +  ?)  >< 

o  o 

o 

g.  (g  =  s(p)  refers  to  electric  vector  polarization  perpendicular 


(iic.x-ftt) 

e  1  ,  where  ft  =  c|k^|  and  is  the  amplitude  in  the  polarization 


(parallel)  to  the  plane  of  incidence.)  Then,  using  Fresnel's  equations 
to  give  the  amplitudes  of  the  reflected  and  transmitted  waves ,  we  write 
the  total  electromagnetic  potential  as  a  sum  of  incident,  reflected, 
and  transmitted  waves  as 


-K  ,  -iftt 


A(x,  t)  =  [A3  (x)  +  (x)  ]  e 


(ID 


where 


■fs  ■> 
A  (x) 


lSelKR( 


;  x>  0 

1+  nra 


o  1  +  nra 


2  iKR  — ibx  v  ^  n 

e  e  n  *  K  ;  x<0 


(12) 


aP(x)  =  - 


tP  iKR 
Ar  e 
o 


i£x  .  nr  a  -i£xN  i?x  nr  a  -i£xNl  ^  n 
n(e  *  + - —  e  )-  cK(e  -~~rr  e  )l  ;  x  >  0 


Kn( 


nr+a 


nr+a 


-fp  2  iKR  -ibx  .  n 

A  - 7 - rr  e  e  (Kn  +bK)  ;  x  <  0 

o  n  (n  +a) 
r  r 


with  a 


cos  et  i  r~i  o  2 

=  - -  “/ K  (1  -1/n2)  +  c  • 


cos  0  C  I" 

A(x,t),  as  defined  in  (11)  and  (12),  is  a  solution  of  Maxwell's  equa¬ 


tions  (8)  . 


•I 


14 


Instead  of  the  potential  A  (x) ,  one  requires,  for  quantization,  a 
set  of  orthonormal  modes.  Define  a  prospective  set  of  modes  as 


n  1  iKR  ,-K 

Vx)  =  e  Vx) 


(13a) 


where 


fj  (R,x)  =NSnx  K- 
KC 


,  i£x  ,  1  nra  -iCx, 

(e  +inr7e  >  ;  x>0 


2  -ibx  _ 

e  ;  x  <  0 


1  +  nra 


(13b) 


fP  (R,x)  =  NP- 


nr  a  _-i£x.  i£x  nr  a  -iCxN 

+v^e  >- ?K(e  );  x>0 


nr(nr+  a) 


(Kn  +bK)  e  lbx 


;  x  <  0 


S  D 

(N  and  lr  are  normalization  constants.) 


.  -v 


The  set  (U.  (x)  }  is  a  suitable  basis  in  which  to  expand  the  electro¬ 

ns 

magnetic  field  resulting  from  an  external  radiation  source  on  the  vacuum 
side.  As  shown  in  Appendix  A,  this  set  obeys  the  orthogonality  rela¬ 
tion 


dx  n2(x)  U+  (x)  U*  (x)  =6  ,  6(lc- kT' )  . 

kg  k*  3*  33 


(14) 


The  modes  (U.  (x)}  are  orthogonal  with  respect  to  the  weighting  factor 
2  kB 

n  (x)  defined  in  (5b)  .  In  Section  (c)  we  find  that  this  form  of  ortho¬ 
gonality  relation  is  crucial  in  determining  the  form  of  the  Hamiltonian 

and  the  commutation  relations  for  the  operator  C+  . 

kS 

We  shall  also  use  in  (c)  the  additional  relation: 


dx  n2(x)U_v  (x) -U^  (x)  =  -  A*  5  (K  +Kf )  6  (  C  -  C ' ) 

kS  k'B'  kg 


(15) 
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This  relation  is  derived  in  Appendix  A,  where  the  coefficient  A->  is 
defined. 


c)  The  quantization  formalism 

To  quantize  the  electromagnetic  field,  we  shall  follow  the 

12 

Lagrangian  formalism,  as  described,  for  example,  in  Haken  .  We  must 
first  find  a  suitable  Lagrangian  that  is  consistent  with  the  desired 
equations  of  motion  and  Hamiltonian.  Since  the  electromagnetic  fields 
(E  and  B)  are  derived  from  the  potential  (A)  ,  we  will  find  the  Lagran¬ 
gian  (L)  to  be  a  functional  of  the  coordinate  A(x, t)  and  its  deriva¬ 


tives.  The  equations  of  motion  are  then 


_d 

dt 


SL 


6L 


SA^Cx,  t)  6Ai(x,  t) 


=  0 


(16) 


where  A_^(x,t)  is  a  Cartesian  component  of  it(x,  t)  and  the  dot  means  8/8t 
Equation  (16)  must  be  consistent  with  the  wave  equation  (8a)  that  was 
previously  derived  from  Maxwell’s  equations.  This  is  so  if  the 
Lagrangian  is  chosen  to  be 


L  = 


8tt 


dx 


3 

'  -  l  [(A.(x'))2  -(VA  (5>)2] 

vZ(x')  i=l 


(17) 


The  constant  of  proportionality  in  the  above  expression  was  fixed 

by  requiring  that  the  Hamiltonian  be  consistent  with  the  classical 

.  12 

field  energy  expression 


-■I 


,  ->2  ,-K  1  ->2  x 

dx  (e  E  (x)  +  —  B  (x) )  . 


(18) 


Using  the  chosen  Lagrangian  (17),  one  may  define  the  canonical 
momentum,  conjugate  to  the  coordinate  A^(x)  as 

6L 


n±(.x,  t)  = 


6Ai(x,  t) 


2  ^ 

=  n  (x)  A^x,  t)  . 


(19) 


. 
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Now,  one  makes  the  quantization  postulate,  namely  that 


[A.  (x' )  ,11 .  (x)  ]  =  iftS^Cx  -x'  ) 
^  J  1 J 


(20) 


where 


**■»  \  P  f -*■»  \  ,19  1 

o..(x-x  )  =  6..6(x-x  )  +  - - - — - - — 

ij  ij  4tt  9x.  9x. 


(21) 


1  J  I  x  ~  x  , 

The  first  term  in  (21)  is  just  the  usual  delta  function,  while  the 
second  term  is  included  so  that  the  commutator  (20)  is  consistent  with 
the  gauge  condition,  =  0. 

The  commutator  (20)  for  the  field  coordinates,  combined  with  the 

expansion  for  ^(x, t)  given  in  (9)  and  the  orthogonality  relations  (14) 

and  (15),  together  imply  the  commutation  relations  for  C->  and  C* 

kB  k' 3 ' 

given  below: 


[VirI  =  Vs(k+4,) 
[Vci?'e-1  =  0=[i’  iB-] 


(22a) 

(22b) 


In  fact,  the  commutators  for  C^.  are  given  by  equations  (22a, b)  only  if 

k3 

A^(x)  and  II (x)  obey  equation  (20). 

One  may  then  express  the  Hamiltonian  in  the  form 


H  =  I 

6 


dk  hn.  (Cl  C+  +  h 

£3  ica  £3  2 


(23) 


where  =c|k|.  To  derive 

iTb 


(23),  substitute  the  expansion  (9)  for  A 


into  the  expression  for  H  given  in  (18).  Then,  using  the  orthogonality 
relations  (14)  and  (15),  and  the  commutators  (22a, b),  the  given  result 
is  obtained. 

Interpreting  CL.  and  CL>  as  the  annihilation  and  creation  operators, 

ka  ka 

respectively,  of  the  photon  states  (k,8)  of  the  electromagnetic  field. 


. 


' 


the  quantization  procedure  is  almost  complete.  Formally, 


one  must 


further  check  that  the  chosen  set  of  modes  {ft*  (x)  }  is  complete.  This 

k6 

is  required  for  the  theory  to  be  self-consistent.  We  have  not,  howeve 
performed  this  calculation. 


d)  Perfect  conducting  solid 

If  the  solid  is  a  perfect  conductor  (infinite  conductivity)  in¬ 


stead  of  a  dielectric  (zero  conductivity),  then  one  may  perform  the 
same  quantization  procedure,  provided  one  uses  instead  the  orthonormal 
modes 


/"N  >Ts  iKR  ,  irx  -irxx  ->  -> 
tU  (x)  =  N  e  (e  s  -  e  ^  )  n*  K 
ks 


x  >  0 


rF  /“N  xtP  iKR 
U .  (x)  =  Ne 

kp 


KS(e1?x+e  1?x)-  d(eltx-  e  iCX)l  ;  x>  0  (24) 


Uv  (x)  =  U  (x)  =  0  ;  X  <0 

ks  kp 


with 


dxlL  (x)U*  (x)  =  6fla, «(£-£') 

kg  k'g'  $S 


dS  tL  (x)%  (x)  =  5  fS(it+£’)<5(c-Cf) 

kg  kg'  SS 


(25) 


The  quantization  done  in  Section  2.2c  then  carries  through  iden¬ 
tically  for  this  set  of  modes. 


§2.3  Hamiltonian  of  the  Gas-Solid-Laser  System 

We  shall  now  construct  the  Hamiltonian  for  the  gas-solid-laser 
system  in  second  quantized  form.  We  assume  that  the  density  of  gas 
molecules  is  small,  both  on  the  solid  surface  and  in  the  gas  phase,  so 
that  gas  molecule-gas  molecule  interactions  can  be  neglected.  The  gas 
molecules  may  then  be  treated  independently. 


. 


§2.3.1  The  static  Hamiltonian 


We  study  one  diatomic  molecule  in  a  vacuum  near  a  semi-infinite 
static  solid  surface,  and  will  ignore  thermal  vibrations  and  the  elec¬ 
tromagnetic  field  for  now.  See  Fig.  5.  Atoms  one  and  two  of  the 
molecule  are  separated  by  the  displacement  t  =  -x  . 


Fig.  5.  Adsorbate  molecule  configuration. 


Now  let  U(t)  be  tbe  interatomic  potential  of  the  molecule,  and 

th 

(x. )  be  the  potential  energy  of  the  i  gas  atom-solid  interaction 
(i  =  1,2).  The  quantum  mechanical  Hamiltonian  for  the  diatomic  gas 
molecule  and  the  static  solid  is 


h  =-  j- 
s  t  2m. 


ft 


1  a5? 


2m 


2  9x: 


^2  +  U(?)  +  V1(x1)  +  V2(x2)  . 


-  2 

Introducing  centre  of  mass  coordinates,  H  becomes 

u  L. 


(27) 


' 


■ 
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ft2  z2  h2  a2 


+  +  V*  +  ^> 


iL 

h 


(28) 


where 


m  =  m1  +  m2  , 

,  —1  -1  -1 
y  =  (m1  +  m2  ) 


x  = 


=  -  (m1S1  +  m252)  , 


■±  ■>  + 

5  =  x2  '  *1  . 

To  simplify  calculations  considerably,  we  choose  to  ignore  the 

interaction  between  atom  2  and  the  solid,  and  thus  set  V0(x+  —  f)  =0, 

2  m£ 

Atom  2  is  sufficiently  distant  from  the  surface  to  render  its  inter¬ 
action  with  the  solid  insignificant.  Then,  defining 

V?)  =  -!^  +  u(l) 


H  (x)  =  +  ) 

l  2m  ^2  1  m^  o 


m 


(29) 


and 


Hres^>  '  V*  -%t>~  V*  -  £  V 

where  £  is  the  mean  interatomic  spacing  of  atoms  1  and  2,  one  has 


st 


1) 

=  H(t) 

+  H 

y 

m 

->•  i 

i  -*  i 

U- 

-g  « 

l«ol  ' 

-k 


res 


(30) 


term  (H  )  in  the  Hamiltonian  will  be  small;  H  is  therefore  treated 
res  res 

as  a  perturbation. 

It  is  desirable  now  to  find  a  diagonal  basis  for  the  operator 

H  (f)  +H  (x) .  H  and  H  can  be  diagonalized  independently,  as 
y  m  y  m 


expressed  by  equations  (31a)  and  (31b). 
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Vt)uv(|)  =  AVUV(^ 

H  (x) $ . (x)  =  e . $ . (x) 
m  J  J  J 


(31a) 

(31b) 


where  u^(0  is  an  eigenvector  of  with  corresponding  eigenvalue  Av 

and  $^(x)  is  an  eigenvector  of  with  corresponding  eigenvalue  £  ^  . 

(For  our  choice  of  V^(x)  and  U(£)  -  Morse  and  harmonic  respectively  - 

the  solutions  to  (31a)  and  (31b)  are  known  in  analytic  form.  See 

Appendix  B.)  Since  H  and  H  are  Hermitian  operators,  the  correspond- 

y  m 

ingly  sets  of  eigenvectors  are  complete  and  may  also  be  orthonormalized. 
Then 


H  (?) 

y 


H  (x) 
m 


<f>.(x)u  (?)  =  (£  .  +A  )$.(x)u  ( I ) 
J  V  J  y  j  V 


(32) 


We  shall  now  adopt  the  formalism  of  second  quantization  to  des¬ 
cribe  the  total  system  of  many  gas  molecules.  The  second  quantized, 
many  particle,  static  Hamiltonian  is  then 


H  = 


dx  d£  1/;  (x,5)H(x,£)tJj(x,£) 


(33) 


We  expand  the  field  operator  in  terms  of  the  eigenvectors  of  (H  +  H^) 


as 


i|j(x,t)  =  l  *.(J)u  (?)  ou 

jfV  3  J 


(34) 


aT(aV  )  is  the  annihilation  (creation)  operator  for  a  particle  in  the 
J  J 

th 

j1"*1  state  of  the  surface  potential  and  in  the  v  molecular  vibrational 
state.  Then,  applying  (32)  through  (34)  and  using  the  orthonormality 


of  {u  }  and  {$.}  leads  to 
n  J 


(H  +  H  )  = 
y  m 


dx  d?  /(x,?)  H  (?)  +  H  (x)  i p(x,t)=  l  E^aT+aT 
Ly  mj  j ,  v  J  J  J 


where  E ;  =  e  .  +  A  »  and 

J  J  v 


(35) 


. 


H. 


res 


dxdt^(x,5)H  (x,f)iKx,5)=  l  Z(j,v,j',v)aV+ay,  (36) 

jj’w'  ^  3 


where 


Z(j,v,jT  ,v’)  = 


dx  df  $*(x)u*(l)H  (x,f)$.f(x)u  , (I) 

j  v  res  j  v 


(37) 


§2.3.2  The  electromagnetic  field  interaction 

In  Section  2.2  the  electromagnetic  field  due  to  laser  radiation 
in  the  absence  of  any  gas  molecules  was  determined,  with  Hamiltonian 


given  by  equation  (23)  as 


H  =  l 

em  £ 


dk  (C+  C_>_  +  — )  . 

kg  kg  kS  2 


(38) 


To  calculate  the  effect  of  the  electromagnetic  field  on  an 
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adsorbed  gas  molecule,  we  use  the  method  of  minimal  coupling  .  The 
kinetic  momentum  p  of  a  point  particle  with  charge  Q  and  mass  m  is 
replaced  by  p  -Q^t,  and 


2 

2m  V 


.,2  2 
-ft  V 


2m  m 


(39) 


Assuming  that  the  molecule  comprises  point  particles  1  and  2  with 
charges  Q  and  -Q,  respectively,  substitution  of  (39)  into  (27)  gives 


H 


st 


H  -  ih  -0-  A(x,  t)  •  — 

st  v  at 


(40) 


We  have  made  the  approximation  that,  since  the  radiation  wavelength  is 
much  larger  than  |£|,  A(x)  changes  very  little  over  £.  Thus 


A(x.)  =  A(x)  +  VA(x).(-iL)  5 
i  hk 

where  i  =  1,2.  The  laser  coupling  term  is  then 
H  =  -  rh  ^  A(x,t)*-4?  , 


which  gives  in  second  quantization 


-  ' 


r 


-K  -K 


dx  df  ^+(.x,f) H^(x,5)ip(x,5) 


=  -iff  ^  I  dx$*(x)l(x,  t)$. ,  (x)» 
y  jj'w'  J  J 

v+  V1 

*  a  .  a. , 
a  j 


df  u*(£)-~  u  ,  (f) 
v  35  v 


(41) 


Now,  assuming  the  adsorbed  molecule  is  localized  at  a  particular  site 
on  the  surface,  we  may  take  A(x, t)  =  1(0, t)  where  the  origin  is  chosen 


to  be  on  the  surface  at  the  adsorption  site.  Then  (41)  becomes 


H 


Z 


=  -Hi  ^  X(.0,t)  •  l  6  ,  fd|  u*(|)  —  u„,(E)«;’o’,  . 


-  •  I  I  JJ 
JJ  WT  JJ  J 


V 


at 


v+  v’ 
i .  a . 

J  J 


(42) 


(This  is  the  so-called  dipole  approximation.)  Expanding  A(x,t)  as  in 
(9)  and  defining 


Y(j,v,j\v')  =-ih  l~-  ^5.., 

*  2eo  ^  JJ  J 


d£  u*(5)-—  U  ,  (5)  , 

v  35  v 


(43) 


one  may  rewrite  the  laser  interaction  term  in  the  Hamiltonian  as 


=  I  Y(j,v,  j',vf)  aT+aT,  l 


JJ  'vv1 


J  J 


dk 


kg 


-iO*  t 
U  +  (0)  e  k6  C+  +h.c. 


kg 


kg 


(44) 


§2.3.3  Dynamics  of  the  solid  -  phonons 

An  adequate  treatment  of  the  solid  must  not  neglect  its  internal 
thermal  vibrations.  We  assume  that  the  only  relevant  effect  of  these 
vibrations  is  to  move  the  surface  adjacent  to  the  adsorbate  molecule 
from  x  =  0  to  x  =  a.  a  is  a  function  of  time  and  may  be  expanded  in 
terms  of  the  normal  modes,  or  phonons,  of  the  solid  as 


a(t)  = 


h 


I 


-+• 

e 


2M  N 
s  s  J 


0) 


b+(t)  +  bj(t) 


(45) 


where  M  =  mass  of  a  solid  molecule,  N  =number  of  molecules  in  solid, 
s  s 


■ 
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J  is  a  multi-index  identifying  the  vibrational  modes  (phonons)  of  the 

solid, each  of  which  has  annihilation  (creation)  operator  b  (b  ),  unit 

J  J 

polarization  vector  e  ,  and  frequency  to  . 

J  J 

To  include  the  effect  of  thermal  vibrations  in  the  Hamiltonian, 
let  the  surface  to  molecule  displacement  now  be  x  -a(t)  in  equation 
(28)  (and  ignore  V^)  to  get 


H  *  H  +  V  (x-a(t)  -  O  ~  V  (x  -  O  . 
st  stl  ml  1  ml 


(46) 


-> ,  „  y  -K 


Expanding  V  (x-a(t) — —  £)  in  a  Taylor  series  to  first  order  in  a(t) 
i  ^1 

gives 


H 


st 


H  -  a(  t)*  -4^  V  (x - —  I)  • 

st  3x  1  m 


(47) 


The  validity  of  this  approximation,  which  includes  one  phonon  processes 
only,  is  discussed  in  references  2  and  15.  We  note  here  only  that  the 
coupling  to  the  phonons  in  the  solid  must  be  weak. 

Defining 


H 


ph 


-  a(t)  V  (x 
3x  1 


r-  O 

mi 


(48) 


and  using  (45),  the  phonon  interaction  term  in  the  second  quantized 
Hamiltonian  becomes 


H  .  = 
ph  J 


-k 


dx  d|"^+(x,f)  Hph^(x,C) 


h 


l 


-> 

e 


J  /v  +  , 


,  L  —  (b‘[(t)  +  b  (t))  l 

2MsNs  J  ^  J  J  jj’wT' 


*  ->•  *  ,-k  3 
dx d£  $  .  (x)u  (O  — 
J  v  3x 


y  -K  V+  V 

X  V  (x  -  £)$.,(x)u  ^) a.  a  , 

1  j  v  J  j  j 


or 


H  -  l 

ph  jj'vv' 


x(j  ,V,  j  '  ,vf  )aY  aj,  •  l  (bl(t)  +  b 


-v 

e 


J  „  f 


J  J'  5  J 


(t)] 


(49) 


24 


where 


X(j,v,jf,v’)  =- 


ft 


2M  N 
s  s 


dx  d£  (x)uv(5)  —  V1(x-  —  O^x^,  (£) 

(50) 


We  have  assumed  here  that  the  phonon  spectrum  of  the  solid  is  essen¬ 
tially  unperturbed  by  the  adsorbed  molecule,  and  hence  aj  commutes  with 
b  .  (This  approximation  excludes  the  possibility  of  heating  the  solid 
via  this  interaction.)  The  Hamiltonian  of  the  solid  is,  in  the  harmonic 
approximation, 

“s -  \ >J  (bI  bj+  i>  •  <51> 

vJ 


§2.3.4  The  complete  Hamiltonian 

We  are  now  prepared  to  construct  the  complete  Hamiltonian  for  the 
dynamical  gas-solid-laser  system.  It  is  presumed  that,  initially,  (for 
t  <0)  equilibrium  conditions  exist.  The  initial  occupations  of  states 

A  A  A  A 

are  assumed  to  be  given  byH+H+H+H  ,  which  is  diagonal  in  the 

y  m  s  em 

basis  we  use.  Then  at  t  =  0  the  laser  is  turned  on,  creating  a  non¬ 
equilibrium  situation.  The  subsequent  time  evolution  of  the  system  is 

A  A  A 

determined  by  the  non-diagonal  interaction  terms  H  ,+H  +H  ,  so  we 
J  o  ph  £  res 

turn  on  these  terms  at  t  =0.^  The  Hamiltonian  may  then  be  written  as 

H  =  H+H+H+H  +  0  ( t)  [H  ,  +  H  +  H  ]  .  (52) 

y  m  s  em  ph  l  res 

Formally,  H  should  be  included  as  a  correction  to  H  for  all 
res  m 

A 

times.  However,  it  will  be  seen  in  the  next  section  that  H  does  not 

contribute  when  the  gas  molecules  are  restricted  to  adsorbed  states. 

This  is  essentially  the  case  for  t  <0,  since  the  occupation  numbers  of 

the  surface  states  initially  follow  a  Boltzmann  distribution  and 

k  T  <<  V  .  (V  is  the  depth  of  the  potential  well  of  the  surface  bond.) 
Boo 


Note  also  that  we  have  as s uae d  that  the  laser  is  kept  on  through¬ 
out  the  complete  desorption  process.  Thus,  when  a  pulsed  laser  is 
used,  our  theory  predicts  a  desorption  rate  that  carries  on  for  only 
the  duration  of  the  pulse. 


§2.4  Calculation  of  the  Transition  Rates 
§2.4.1  Application  of  the  perturbation  theory 

In  Section  2.1  we  proposed  to  describe  the  desorption  process  by 
the  master  equation 


nj(t)  = 


vv’  vT 


l  R::,n'(t) 

( j '  ,vf)  £(j,v)  32 


rv:v+  v  rv,.v 

3  3  v'£v  CJ 


nV(t)  . 
3 

(2) 


V1 V  y’ y 

We  will  calculate  the  contributions  to  the  rates  R. . .  and  R  .  due 

3  3  cj 

A  A  A. 

to  the  interaction  terms  H  ,  H  ,  ,  and  H  individually.  Using  Fermi’s 

l  ph  res  •  ® 

Golden  Rule,  (second  order  tine -depen dent  perturbation  theory)  the 
single  particle  transition  rate  from  state  a>  to  b>  is  given  by-' 


-  =  4r.  <b  H 

D-<-a  h  m 


Lnt  fi(Eb  -  Ea> 


(53) 


where  H.  is  the  interaction  term  in  the  Hamiltonian  and  E,  ,  ,  is  the 
mt  (a,b) 

energy  level  of  state  (a,b).  Equation  (5  3)  is  valid  in  the  long-time 

limit,  and  provided  that  the  interaction  (H.  )  is  sufficiently  weak. 

in  t 

The  state  of  a  single  adsorbed  gas  molecule  is  described  by  its 
surface  bond  state  (j)  and  internal  vibrational  state  (v)  which  have, 
in  our  approximate  model,  a  separable  state  function  j>  v>.  The  state 
of  the  solid  is  characterized  by  its  phonon  distribution  {nr'n},  pre- 

vJ 


viously  assumed  to  be  basically  unperturbed  by  the  gas  molecules  and 
laser.  The  laser  electromagnetic  fields  are  described  by  the  photon 


1 1*1:;.  liiC  -i  £  .  S 


26 


distribution  {n^m} ,  independently  of  both  phonons  and  gas  molecules. 

The  complete  one-gas-molecule+solid+laser  system  may  then  be  described 
by  the  separable  state  function  |  j >  |  v>  [  {nP^}>  |  {n^.™}> .  We  can  now  apply 
(53)  to  our  system  to  calculate  the  single  particle  transition  rates. 
According  to  (53),  the  transition  rate  from  (j',v?)  to  (j,v)  due 

a 

to  the  laser  term  (H  )  is 

Xy 


w’  _  2ir 
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x  6  E 


k3  '  * 

-  E 


j  ,  v,  (nPh> ,  {n^}  j  ’  ,  v’  ,  {nPM '  ,  {nem}  ’ 


(54) 


Substituting  expression  (44)  for  H  in  (54)  and  noting  that  H  is  inde- 

X/  36 

pendent  of  phonons,  we  get,  after  a  few  simple  algebraic  steps. 


w'  _  2jt  r- 

j  j f  “  ft  \ 
jj  3 


dic?rr  |Y(j,v,  j' ,v')  (0)  |2 

kg 


'nt 


n^.m  6(EV  -EV[-  fiSK)  +  (nJm+l)6(EV-  EV!  +  hO*) 
L  k3  J  J  k  k3  J  J  k 


(55) 


Similarly,  we  substitute  (49)  for  H  ,  into  (53)  to  get  the  phonon  assis- 

ph 

ted  transition  rate. 


„wT  2tt  v  1  it/ .  . »  .  N  1 2 

Pjj,-T  ^  X(3’V’J  ’v  ),ejl 

vJ 

X  [nPh6(EV-  Ev'  --hoiT)+  (nph  +  l)  (EV-  EV|  +liai  ) 
J  1  l  J  J  11  J 


(56) 


(We  may  replace  the  bound  state  index  j  with  q  to  get  the  transition 
rate  to  a  continuum  state.)  Finally,  the  transition  rate  due  to  tunnel- 

A 

ing  via  H  is  calculated  by  substituting  (36)  into  (53)  to  get 

1T6S 


q™:  =f  izcq.v.j-.v')!2 


(57) 


We  assume  that  the  energy  levels  E^  are  non-degenerate  for  all  bound 


states.  This  is  plausible  since  e.  and  A  take  on  discrete  and  un- 

J  v 
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related  values,  so  it  is  highly  improbable  that  eY  -Evf  =0  if  both  j* 

t  ^  1 

and  j  are  bound  states.  Therefore  qYY,  =  0.  If  the  final  state  is  in 

33 

the  continuum,  it  is  then  possible  that  EV-eY,  =0,  and  hence  QVV,  may 

q  3  qj 

be  non-zero.  It  should  be  remarked  that  Lucas  and  Ewing^  have  consi¬ 
dered  desorption  due  to  this  tunneling  term  alone.  We  will  later  see 
how  their  theory  fits  into  ours. 

The  calculation  of  the  transition  rates  LYYr  ,  pYY,  ,  and  QW, 

JJ  JJ  qj 

involves  the  evaluation  of  the  matrix  elements  Y( j , v, j ’ ,  vf ) ,  X(j,v,j',v') 


and  Z( j, v, j ' , v' ) ,  respectively.  To  simplify  the  calculations  in  the 


following  sections,  we  assume  that  the  axis  of  the  adsorbed  diatomic 

molecule  is  perpendicular  to  the  surface.  We  furthermore  restrict  atoms 

1  and  2  to  motion  along  this  axis,  and  take  the  wave  function  for 

the  molecule  to  be  one  dimensional,  hence  $ . (x)  =  $.(x).  This  one- 

3  3 

dimensional  approximation  is  really  more  appropriate  if  the  molecule  is 
free  to  move  on  the  surface,  rather  than  highly  localized  as  we  assume 
here.  We  do  not,  however,  expect  this  approximation  to  affect  the 
results  significantly  in  a  qualitative  way^ . 


§2.4.2  Y(j,v,j',v')  and  lYY, 

2  J _ 

To  calculate  L^TY,  ,  we  first  evaluate  |  Y(j ,  v,  j  ’  ,vf )  *U->  (0)|^  and 

J  J  k3 

nY.m  ,  and  substitute  these  expressions  into  equation  (55)  for  lYY,  • 
kg  33 

Y( j , v, j f , v' ) ,  given  in  (43),  is  now  rewritten  in  the  one  dimensional 
approximation  (£=£ n,  x=  xn)  using  bra-ket  notation  as 


-> 


Y(j,v,jf,v')  =-ifi  ^  6  , 

<ze0  q  33  ) 
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=  -ifi 
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- —  x  <  v  —  v  >n  6  .  .  ,  . 

2e  ^  H  JJ 


(58) 
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Since  the  molecular  vibrations  are  assumed  to  be  simple  harmonic, 
we  make  use  of  the  annihilation  (creation)  operators  a  (a  )  of  harmonic 
vibrational  modes  to  express 


yoa 


<v 


v'>  = 


2h 


V  I  ti  f 

< v  a  -a  v  > 


yoj  v 


2h 


Ar  +1  <5  ,  ,  —  /v  6  f 
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(59) 


where  co^  is  the  fundamental  molecular  vibrational  frequency.  Combin¬ 
ing  (58)  and  (59)  gives 

|Y(j,v,j’,v’)\(0)|2  = 


(v+1)  5  ,  +v  5  , 

v  , v+1  v  ,v-l 


(60) 


yQ  y 

where  U->.  is  IL>  (x)  evaluated  as  x^-0  with  x  >  0. 

kg  ke 

(Note  that  only  p-polarized  photon  modes  contribute  to  (60).) 
Substituting  (60)  into  (55)  and  rewriting  the  delta  functions  slightly 


gives 
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(61) 


It  is  convenient  to  introduce  polar  coordinates  in  (61)  as  k  -*kk  to  get 


_  q  ir  5 
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(62) 


Suppose  the  laser  radiation  is  strongly  concentrated  about  one 


^  6m 

direction,  k  ,  so  n  ,  is  sharply  peaked.  The  slowly  varying  term 
°  kkp 

|n*U°  |  may  then,  to  a  good  approximation,  be  removed  from  under  the 


kfcp 
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integral  over  solid  angle  k.  The  integral  over  k  is  easily  evaluated 


using  the  delta  function,  leading  to 
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em 


It  is  useful  to  express  n^  in  terms  of  an  easily  measured  quantity, 
such  as  the  incident  laser  intensity  (I)  which  we  take  to  be  the  ampli¬ 
tude  of  the  incident  Poynting  vector.  For  plane  weaves  in  a  vacuum  one 
has 


I  = 


,  -v 
S 


ave 


o  | -vi  .-v  .  ,2 

=  hr  lE  (x,t) 


ave 


(64) 


where  ave  means  the  time  average  and  superscript  "i"  means  that  only 
the  incident  component  of  the  wave  is  included.  Considering  now  the 
contribution  due  to  the  incident  portion  of  a  single  mode  (k6), 
equations  (7)  and  (9)  in  (64)  yield 


X->  =  c|u£  \2[m+  (njm  +  \)  ] 

ke  ks  k6  k6  2 


(65) 


Note  that  this  is  just  the  energy  density  times  the  wave  velocity, 
Ignoring  the  zero  point  energy  contribution,  (65)  may  be  rewritten  as 
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em 
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(66) 


em 


(This  approximation  is  valid  since  we  find  Uj-v  >>  1.) 

A  .  | 

Since  I .  is  sharply  peaked  about  k  and  IL. 

Ic6  °  kg1 

varying,  we  get  approximately 
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(67) 


To  arrive  at  the  second  equality  above,  we  have  taken  the  total 
laser  intensity  to  be  I  with  normalized  frequency  distribution  g(fi) . 

The  factor  of  1/2  is  due  to  the  fact  that  only  p-modes  contribute  to  the 
interaction,  but  we  assume  that  the  laser  radiation  is  unpolarized. 
Substituting  (67)  into  (63)  we  get 
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r  photon  ^  ('Stimulated') 
^absorption'  ^  emission  ' 


(spontaneous  emission)  (68) 


To  proceed  further,  details  of  the  laser  radiation  must  be  speci¬ 
fied.  At  normal  laser  intensities,  the  last  term  in  (68)  -  due  to  spon¬ 
taneous  photon  emission-  can  be  safely  ignored.  We  shall  also  take  the 
frequency  distribution  to  be  Lorentzian,  centred  on  the  molecular  vibra¬ 
tional  frequency,  co  .  That  is, 

g(S2>  =  7  - - 2  (69) 

(a  -  wv)  +  (y/2 r 

where  y  is  the  half-width  of  the  frequency  distribution.  If  the  laser 
is  in  continuous  wave  operation,  then  y  is  the  line  width  of  the  radia¬ 
tion  associated  with  decay  processes  of  time  constant  t^.  On  the  other 
hand,  a  laser  in  pulsed  operation  -  "Q-switched"  with  pulse  duration 
t  ~  100  ns  -  will  have  frequency  spread  given  by  the  Fourier  transform 


'  i*. 
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of  the  pulse  shape,  and  hence  y  ~t  \  (We  assume  that  the  frequency  spread 
due  to  the  pulse  shape  dominates  over  all  other  sources  of  line  broad¬ 
ening;  that  is,  t  <<x  .)  In  either  case,  the  frequency  distribution 

P  ~ 

may  not  be  Lorentzian,  but  the  error  introduced  here  is  small.  Equa¬ 
tion  (68)  then  becomes 
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For  a  pulsed  laser,  we  let  the  pulse  fluence  be  F„  =  I  t.  (In 

J6  O 

fact,  in  the  case  of  continuous  wave  laser  operation,  we  may  use  the 
term  "fluence”  as  defined  above  to  mean  the  intensity  times  inverse 
line  width. ) 

For  our  dielectric  solid,  equations  (13a, b)  give 
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(71) 


(Note  that,  as  n^  1,  (71)  sin  9.  This  is  just  what  one  expects  in 
an  infinite  vacuum  without  the  solid.)  The  final  expression  for  the 


transition  rate  is 
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From  the  form  of  L^,  above,  it  is  clear  that  the  photon  interaction 

JJ 

induces  one— step  transitions  between  vibrational  (v)  states,  but  does 


no 


t  affect  the  surface  state  (j). 


As  an  aside,  we  note  that,  if  the  solid  were  instead  a  perfect 
conductor,  one  must  use  equation  (24)  to  replace  the  form  factor  (71) 
in  (72)  with 


2 


(73) 


Both  form  factors  give  zero  contribution  at  perpendicular  incidence 
(9  =0).  While  the  form  factor  for  the  perfect  conductor  is  largest  at 
grazing  incidence,  the  dielectric  form  factor  is  maximized  at  an  inter¬ 
mediate  angle.  See  Fig.  6  in  which  (71)  and  (73)  are  compared  graphi¬ 
cal  ly . 
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Fig.  6.  Angular  dependence  of  the  vibrational 
transition  rate. 
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§2.4.3  The  phonon  interaction;  X  and  P, 
i)  Bound  state-bound  state  transitions 

We  shall  now  calculate  the  transition  rates  due  to  the  phonon 

A 

interaction,  H  First,  consider  transitions  among  bound  surface 
states  only. 

f 

The  phonon-assisted  transition  rate,  Peof,  given  in  (56),  depends 

— ^ 

on  the  matrix  element  X( j , v, j ’ ,  v* )  .  In  our  one-dimensional  theory  (50) 
becomes 
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We  choose  to  use  a  Morse  surface  potential  extending  perpendicularly 
to  the  surface .  Hence 
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). 


(75) 
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This  potential  has  a  maximum  depth  of  at  x^  and  a  range  given  by  y 
Provided  we  have  xq  >  y  so  increases  rapidly  for  x  <  0,  we  expect 
this  potential  to  be  realistic^".  In  Appendix  B  the  energy  eigenvalues 
and  wave  functions  of  this  potential  are  listed. 

h 


Expanding  (x - C)  in  a  Taylor  series  about  x  in  (74) 
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(76b) 


and  cf>.(£-£  )  =  y  2$.(x).  These  bound  state  Morse  wave  functions  are 
5  ^  J 

given  in  Appendix  B,  and  I  . 

£  J  J 


is  evaluated  in  Appendix  C.  The 
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dimensionless  integration  variable  in  (76b)  is  4  =y(x-x  — —  £  ).  We 

o  m1  o 

choose  the  lower  limit  of  integration  -  equivalent  to  x=0  -  to  imply 
that  the  gas  particles  do  not  penetrate  the  solid.  Since  the  v  states 

are  those  of  a  harmonic  oscillator,  we  express  the  coordinate  £  -  £q 
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2 — ]  <<  1,  so  we  may  write,  to  a  good  approximation. 
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where  v  (v  )  is  the  greater  (lesser)  of  (v,v’). 


To  perform  the  sum  in  (56),  we  must  know  the  phonon  distribution 
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at  the  solid  surface.  This  has  been  calculated  by  Ezawa  ,  and  used  in 

17 

previous  desorption  calculations  by  Goldys,  Gortel  and  Kreuzer 
However,  they  found  that  the  numerical  results  differ  only  slightly 
from  the  parallel  calculation  done  using  the  bulk  Debye  phonons  appro¬ 
priate  for  an  infinite  solid.  We  will  simply  use  bulk  phonons  here  and 
expect  our  results  to  be  qualitatively  correct.  The  phonon  index  (J) 
then  describes  the  phonon  momentum  (p)  and  polarization  (a) .  For  the 

Debye  model,  we  may  substitute 
3N 


^D 


l 

J 


03 


03  d03  , 
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which  absorbs  the  n.e  factor  in  (56). 19  Using  (78),  and  splitting 
the  transition  rate  into  parts  involving  phonon  emission  and  absorp¬ 
tion,  we  get 


w*  vvf  w* 

P. =  p;;f  +  ,  p_t 

JJ  em  jj  abs  jj 


(79) 


where 
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(81) 


The  6-functions  ensure  that  the  initial  and  final  molecular  states  (j,v) 

differ  by  at  most  one  Debye  energy.  The  phonons  obey  Bose  statistics, 
ph 

so  n  is  the  Boson  occupation  number  in  thermal  equilibrium. 

•  It  is  useful  to  scale  all  rates  by  the  Debye  frequency,  to^,  and 
introduce  the  dimensionless  quantities. 
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and 
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TiY' 


(82) 


(a  is  approximately  equal  to  the  number  of  bound  states  in  the  surface 
potential.  See  Appendix  B.)  The  energies  are  all  scaled  by  the  Debye 
energy,  -tfo>  . 

Using  the  results  of  Appendix  C,  (80)  and  (81)  may  be  rewritten  as 
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where 
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0). 
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and  6  =—  -  fundamental  molecular  vibrational  frequency  in  units  of  co  : 

v  wD  H  y  D 


j>(j<)  is  the  greater  (les  ser) of  (j,  j  ’ )  ; 


£j'  ~  £j  =  v)  6V  +  7  (i'~3H2cr0~j  -  j'-i); 
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Notice  that  the  transition  rates  decrease  rapidly  as  |v-v' |  increases 

because  of  the  small  factor  ( (4m2/m^) /r6 )  <<  1.  The  collection  of 

->•  20 
gamma  functions  from  X(j,v,j',vf)  involving  j  and  j'  has  been  shown 

to  strongly  favour  transitions  where  j  and  jT  differ  by  one.  This 

encourages  a  cascade  of  step  by  step  transitions  through  the  bound 


s  tates . 


ii)  Bound  state  to  continuum  transitions 

The  total  phonon- assis ted  transition  rate  from  a  particular  bound 
state  (j',v*)  to  any  continuum  state  with  vibrational  index  v  is 
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I 
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qj 


(85) 


where  the  surface  state  label  q  is  the  (one-dimensional)  wave  vector  of 
the  free  (desorbed)  particle.  Using  (56)  we  have 
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Proceeding  similarly  to  above  in  (i) ,  we  use  bulk  Debye  phonons 
and  pass  to  the  continuum  in  q  to  get,  after  some  computation. 
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The  integration  variable  (x)  in  (88)  and  (89)  is  the  kinetic 
energy  of  the  desorbing  particle,  scaled  by  the  Debye  energy  (fi^p) . 
The  limits  of  integration  combined  with  the  theta  function  allow  de¬ 
sorption  to  a  state  (q,v)  via  phonon  absorption  if  the  available 
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energy  (e  ,-vfi^)  is  with  1  (Debye  energy)  of  the  continuum  (e=0).(That 
v» 

is,  e.,-v6  +1  >0.)  Desorption  with  phonon  emission  can  occur  if 

v  * 

£j»-v6v>0.  The  form  of  the  factors  in  (87)  to  (89)  is  similar  to 
that  in  (83)  and  (84)  for  bound  state  to  bound  state  transitions.  The 
matrix  element  X(q,v,j',vf)  encourages  transitions  in  which  the  v  index 
varies  by  as  little  as  possible  and  in  which  j  is  as  close  to  the  con¬ 


tinuum  as  possible. 


§2.4.4  The  tunneling  transition; 

■  . .  — — , — — . — ^  J 

As  explained  in  Section  2.4.1,  the  tunneling  process  due  to  H 

1T6  S 

can  cause  bound  state  to  continuum  transitions  only.  In  analogy  to  the 
phonon  assisted  rates  of  the  last  section,  we  let  the  tunneling  rate  to 
the  continuum  be 
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(90) 


(To  arrive  at  the  second  equality,  we  have  substituted  equation  (57) 
for  0W;.) 

■qj 

In  one  dimension,  and  using  a  Taylor  expansion  of  H  (29  )  in 


res 


Z(q,v,j',v')  (37)  leads  to 

2V  «  l 
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(91) 


This  equation  is  similar  in  form  to  (76a)  for  X( j , v, j ’ , v' ) .  The  inte¬ 
gral  I  .  is  the  same  as  I . .  . ,  in  (76b),  except  that  the  bound  state 
wave  function  for  j  is  replaced  by  the  continuum  state  function  for  q. 
This  integral  is  very  similar  to  one  calculated  for  reference  1,  so  only 
the  result  is  quoted  in  Appendix  C. 
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Performing  on  (91)  the  same  manipulations  that  took  (76a)  into 
(78)  in  Section  2.4.3  and  substituting  into  (90)  after  evaluating  the 
delta  function  gives 
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(92) 


Although  we  express  the  rate  in  (92)  in  terms  of  the  parameters  go 

and  r  for  convenience,  it  must  be  remembered  that  no  phonon  processes 

are  involved  here.  Rather,  this  process  involves  a  direct  tunneling 

of  energy  between  vibrational  (v)  states  and  surface  states  (j).  The 

remarkable  similarity  of  the  components  of  (92)  and  P^T,  (87-89) 

c  j 

is  due  to  the  very  similar  forms  of  the  matrix  elements  Z  and  J.  (^ 
involves  the  5,^  derivative  of  the  surface  potential  while  Z  contains 
the  (£-l)St.) 

This  completes  our  use  of  the  microscopic,  quantum-mechanical 
theory.  We  can  now  construct  the  full  transition  rates  that  appear  as 
coefficients  in  the  matrix  of  rate  equations  in  (2)  or  (4) .  The  bound 
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state  to  bound  state  transition  rate  is 
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JJ  JJ  JJ 


(93) 


while  the  bound  state  to  continuum  rate  is 


w *  w*  w* 
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§2.5  Solution  of  the  Master  Equation 


To  study  the  time  evolution  of  the  adsorbate,  we  will  solve  the 


rate  equations 


-y.  .  -B-  -y.  . 

n(t)  =  R  n(t)  , 
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subject  to  the  initial  conditions 


— ( E Y—  y  )/k  T 
n^(0)  =  N(0)  e  J  g  B 


N(0)  is  the  initial  population  of  the  adsorbate  and  y  is  the  chemical 

g 

potential  of  the  gas  .  (The  initial  distribution  of  occupation  numbers 
in  thermal  equilibrium  is  given,  for  low  particle  densities,  by  the 
Boltzmann  distribution  above.) 

We  follow  the  method  of  reference  1  and  solve  (4)  by  diagonaliza- 

< "  >  -4->-  ->i 

tion  of  the  matrix  R.  Suppose  R  has  eigenvalues  A_^  and  eigenvectors  e  , 
It  is  then  shown  in  reference  1,  by  basic  linear  algebra,  that  the 


adsorbate  population  is  given  by 
m/*-\  nY(0)  A.t 
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where 
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Now,  let  the  eigenvalue  of  smallest  magnitude  be  X  .  Then,  as  t  -+  °°  y 


N(t)  V 
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where  X  <0.  That  is,  at  large  times,  the  desorption  rate  may  be  char¬ 


acterized  by  the  smallest  eigenvalue,  A^,  if  s^  is  not  too  small.  In 
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practice,  we  find  by  numerical  calculation  that  is  usually  very 
close  to  unity  and  all  other  s^'s  are  normally  much  smaller  in  magni¬ 
tude.  In  this  case  the  relaxation- time  form  given  in  (96)  describes 
the  desorption  process  well.  Often  we  find  a  few  other  s^’s  are  signi¬ 
ficant  too,  in  which  case  the  desorption  process  involves  several 
different  time  scales. 

The  determination  of  the  s^  and  A_^  values  is  a  lengthy  numerical 

calculation,  involving  the  diagonalization  of  the  rather  large  matrix 

R,  and  the  inversion  of  the  matrix  of  eigenvectors.  To  keep  the  matrix 

size  finite,  the  spectrum  of  vibrational  levels  is  cut  off  at  an  upper 

limit  (v  ).  Vibrational  states  are  then  restricted  to  v  =0,1,2,..., 
max 

v  .  Such  a  cut-off  is  not  in  principle  unreasonable.  The  inter¬ 
max 

atomic  potential  of  a  real  molecule  is  in  fact  anharmonic,  so  the 
higher  excited  states  are  not  evenly  spaced.  A  laser  tuned  into  the 
evenly  spaced  low  energy  levels  will  therefore  be  unable  to  resonantly 
excite  the  higher  levels.  The  exact  level  at  which  this  cut-off  occurs 
will  depend  on  the  laser  line  width  and  on  the  degree  of  anharmonicity. 

In  the  next  section,  we  will  perform  the  numerical  calculation 
outlined  above  on  several  systems  of  specific  interest. 
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§3.  NUMERICAL  RESULTS  AND  DISCUSSION 


§3.1  Characterization  of  our  Photo desorption  Model 

The  behaviour  of  our  photodesorption  model  must  now  be  examined, 
before  theoretical  predictions  may  be  compared  with  experimental  re¬ 
sults  as  detailed  in  the  next  section.  Our  theory  is  first  applied  to 


a  se 


t  of  model  gas-solid  systems  for  detailed  study.  The  system  para¬ 


meters  are 


i)  the  fundamental  vibrational  frequency  (oo  )  of  the  adsorbate 


v 


molecule,  its  atomic  masses  (m^,!!^)  and  dipole  charge  separa¬ 


tion  (Q), 


-1. 


ii)  the  depth  (Vq)  and  range  (y  )  of  the  surface  potential, 

iii)  the  dielectric  constant  (n  ),  molecular  mass  (M  )  and  Debye 

IT  S 

temperature  (T  )  of  the  solid,  and 

iv)  the  laser  "fluence"  (F.  =  intensity  x  inverse  line  width)  ,  and 


incidence  angle  (9  ). 

(Note  that  we  will  replace  V  and  T  with  the  scaled  parameters,  a  and 

o  D  o 

r  as  defined  in  equation  (82).) 

The  system  parameters  are  chosen  to  model  the  physisorp tion  of 

CH^F  on  an  NaCl  surface,  insofar  as  the  experimental  values  are  known. 

The  CH^F  molecule  is  believed  to  adsorb  with  the  F  atom  closest  to  the 

surface  and  the  CH^  complex  sticking  outwards.  (The  laser  excites  the 

C-F  bond.)  The  adsorbate  molecule  and  solid  parameters  (i  and  iii 

above)  are,  with  the  exception  of  Q,  generally  all  known.  The  value  of 

Q  is  estimated  from  free  space  measurements  of  the  C-F  bond  dipole 

21 

moment  in  the  gaseous  CH^F  molecule  .  While  it  is  known  how  the  C-F 
bond  vibrational  frequency  changes  as  the  molecule  is  brought  up  against 
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the  surface,  the  effect  on  0  is  unmeasured.  The  error  in  Q  introduced 
here,  however  large,  will  affect  our  results  by  simply  rescaling  the 
fluence  axis  of  the  graphs. 

There  is  also  some  difficulty  in  describing  the  surface  potential. 

The  well  depth  (V  )  may  be  readily  estimated  from  the  heat  of  adsorp- 

o 

3  -1 

tion  as  determined  by  thermal  kinetics  experiments  .  The  range  (y  ) 

of  the  potential  is  not  so  easily  determined  from  experimental  data. 

We  therefore  choose  three  different  values  for  y  and  investigate  the 

desorption  kinetics  in  each  case. 

In  Table  1,  parameters  for  the  three  systems  are  listed.  We  let 

-1  ° 

Vq  =  22.5  kJ/mole  and  y  =.35,  .25,  and  .15  A  in  systems  A,  B  and  C 

respectively.  While  the  lower  two  y  ^  values  are  somewhat  smaller  than 

Table  1.  System  parameters  for  CH^F  on  NaCl . 


System 

-1  ° 
y  (A) 

ao 

r 

A 

.35 

22.8 

53.9 

B 

.25 

16.3 

27.5 

C 

.  15 

8.15 

6.9 

Other  parameters  are: 


T  =  281  K  , 
D 


9 


n  =  1.55  ,  6  =  4.98, 

r  v 

0  =  1  el.  charge. 


e 


60° 
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what  one  might  expect  to  find  physically,  these  three  systems  show  the 
varied  results  attainable  using  our  model.  The  energy  spectra  of  the 
corresponding  surface  potentials  are  shown  in  Fig.  7. 

By  way  of  introduction,  we  first  consider  system  A,  which  is 

-1  o 

characterized  by  the  parameters  cr^  =22.8,  r  =53.9,  and  y  =.35  A. 

The  energy  spacing  between  adjacent  surface  states  is  always  less  than 
one  Debye  energy,  so  purely  thermal  desorption  via  one-phonon  processes 
is  possible. 

Using  the  method  described  in  Section  2,  the  desorption  rate  is 

* 

calculated  as  a  function  of  fluence  and  temperature  .  Then  log  rate  is 
graphed  against  log  fluence  for  various  temperatures  in  Fig.  8;  the 
data  is  redisplayed  in  log  rate  versus  inverse  temperature  graphs  in 
Fig.  9.  On  the  latter  graph  we  include  an  estimate  for  purely  thermal 
desorption,  without  the  laser  or  the  vibrational  coupling.  This  is 
accomplished  by  including  only  the  phonon-assisted  rates  an<^  ^cj* 

in  the  calculation,  which  is  precisely  the  thermal  desorption  calcula¬ 
tion  of  reference  1.  The  thermal  desorption  graph  is  a  straight  line, 
indicating  an  Arrhenius  form  for  the  rate  (see  Introduction) ,  with  the 
heat  of  desorption  approximately  equal  to  the  depth  of  the  lowest  bound 
state . 

The  photo desorption  graphs  show  enhanced  desorption  rates  at  low 
temperatures.  At  100  K,  photo desorption  occurs  at  a  rate  100  times 


*We  set  vmnv  =2  for  the  calculations  in  this  section.  A  molecule  in  the 

surface  ground  state  (j  =0)  needs  a  minimum  vibrational  excitation  of 

v=2  to  bring  its  total  energy  into  a  degeneracy  with  the  continuum 

states.  We  find  that  the  photodesorption  rates  are  greatly  hindered  if 

vmax  -^s  reduced  below  this  level.  Increasing  vmax  above  2  results  in  a 
slight  enhancement  of  the  rates,  with  no  qualitative  difference.  We 
prefer  to  keep * * * *  vmax  = 2  here  in  order  to  reduce  the  complexity  of  the 
numerics  to  a  minimum. 


,  -ferns  b©  bfi  •  eiasfc-jsotTr  :r.q  &/13  %£tio  gnJ  nrls/rf.  b&ifaUqmooy* 
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Fig.  7.  Energy  levels  for  the  CH^F/NaCl  systems  A,  B,  and  C. 
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Fig.  8.  Log  rate  versus  log  fluence  for 

system  A  at  various  temperatures. 
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Log  rate  versus  inverse  temperature  for  system  A  for 
different  fluence  values.  (Fluence  is  in  units  J/m  .) 
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larger  than  the  thermal  counterpart,  provided  the  laser  fluence  is 

-2  2 

greater  than  10  J/m  .  (See  Fig.  9).  At  50  K  the  laser  enhancement  is 

of  many  orders  of  magnitude.  The  fluence  dependence  saturates  for 
2 

F^  >  1  J/m  .  This  is  also  seen  in  Fig.  8  where  rate  versus  fluence 
curves  level  off  to  the  right.  Note  that  the  T=200,  250K  graphs  show 
little  fluence  dependence;  thermal  desorption  dominates  in  this  regime. 
(The  Debye  temperature  of  NaCl  is  1^=281  K. )  We  shall  now  move  on  to 
study  system  B,  for  which  a  more  detailed  analysis  is  desirable. 

As  shown  in  Fig.  7,  system  B  has  fewer  bound  states  since  the 
range  of  the  surface  potential  has  been  lessened.  As  a  result,  the 
bottom  three  surface  states  are  separated  by  greater  than  one  Debye 
energy;  thermal  desorption  alone  cannot  occur  from  the  ground  state. 

The  desorption  data  given  in  Figs.  10  and  11  for  this  system  is  rather 
complex,  so  in  the  next  seven  paragraphs  we  investigate  the  fluence  and 
temperature  dependences  with  reference  to  our  microscopic  model. 

The  log  rate  versus  inverse  temperature  graphs  in  Fig.  11  are 
qualitatively  similar  to  those  for  system  A,  except  that  the  high  tem¬ 
perature  rates  roll  off  just  before  meeting  one  another  on  a  common 
"thermal  desorption"  line.  This  effect  will  be  discussed  shortly. 

Before  embarking  on  the  main  analysis,  it  is  instructive  to  modify 
the  rate  equations  temporarily  by  deleting  all  bound  surface  state  tran¬ 
sitions  j->j';  that  is,  PVV,  =0.  The  resultingly  small  desorption 

J  J 

rates  are  graphed  as  dashed  lines  in  Fig.  11.  In  this  very  restricted 

case,  the  only  desorption  path  is  via  (j,v)  =(0,0)  ->-(0,1)  ->(0,2)  ->(c,0). 

The  last  step  in  the  path,  (0,2)  ->  (c,0) ,  is  the  rate  determining  step 

for  fluence  in  the  range  F^  >  10  J/m  .  The  transition  rate  for  this 

step,  which  is  Q°l  +P°^,  turns  out  to  be  only  slightly  temperature 
r  *  ^c0  cO 
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Fig.  10.  Log  rate  versus  log  fluence  for  system  B  at  various 
temperatures . 


Fig.  11.  Log  rate  versus  inverse  temperature  for  system^B  at 
different  fluences.  (Units  of  fluence  are  J/m  .) 
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dependent.  This  is  because  is  temperature  independent,  while  the 
phonon  density  of  states  in  P^q  depends  on  temperature;  see  equations 
(87) — ( 89) .  The  transition  rate  from  (0,2)  to  (c,0)  is  thus  independent 
of  temperature  at  low  T,  but  increases  slightly  as  the  Debye  tempera¬ 
ture  is  approached.  The  corresponding  desorption  rates  shown  in  Fig. 

-3  2 

11  have  just  this  form  for  fluences  F  >  10  j/m  .  However,  for 

X/ 

much  smaller  fluence  the  step  from  (0,v)  to  (0,v+l)  is  slower  and 

_6  2 

becomes  rate  determining.  For  example,  at  F^  = 10  J/m  ,  the  desorp- 

-4  -1 

tion  rate  in  Fig.  11  is  =10  s  at  all  temperatures.  This  is  equal 

to  the  v  =  0  to  v  = 1  transition  rate  L^?.  We  shall  now  turn  our  atten- 

33 

tion  back  to  the  full  system  of  rate  equations, 
w* 

With  P..,  included,  much  more  complex  behaviour  is  observed  in 
JJ 

Fig.  11.  To  assist  in  our  understanding  of  the  underlying  processes, 
we  display  the  numerical  values  of  just  a  few  typical  transition  rates 
in  Fig.  12.  The  vertical  steps  at  the  right  side  of  this  diagram 
correspond  to  the  desorption  path  outlined  in  the  last  paragraph.  It 
is  clear  from  Fig.  11  that  this  path  does  not  normally  contribute  sig¬ 
nificantly  to  desorption  since  the  full-system  rates  are  mostly  much 
greater  than  those  (dashed  lines)  predicted  from  this  single  path  alone 
From  Fig.  12  it  is  evident  that  vib rationally  excited  molecules 
can  transfer  -  via  a  phonon  interaction  -  to  excited  surface  states  by, 
for  example,  (j,v)  =(0,1)  -►(5,0).  The  molecule  can  then  step  up  and 
down  among  the  surface  states  very  rapidly  with  phonon  assisted  transi¬ 
tion  rates  pYY, .  This  "thermal  cascade"  is  included  at  the  left  side 
JJ 

of  Fig.  12.  If  the  molecule  energy  (ej)  is  positive,  desorption  may 
occur  by  a  transition  to  the  continuum.  (See  the  top  of  Fig.  12).  The 
transition  rates  into  the  continuum  from  a  highly  excited  bound  state  - 
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Fig.  12.  Typical  transitions  occurring  in  system  B.  The  columns 
of  bound  states  (j)  correspond  to  vibrational  levels 
v  =0,1,2.  (Rates  are  given  for  T  =  50  K,  in  units  s--*-.) 


■ 
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j  =15,  v-0  or  even  j  =5,  v=l  for  example  -  are  much  faster  than  from 

the  surface  ground  state,  j  =0,  v=2.  This  is  because  the  rates  PV.V 

cj 

V*  V 

and  Q  .  favour  transitions  between  close  surface  states  and  adiacent 

cj  J 

vibrational  states,  as  explained  in  Section  2.  The  actual  desorption 
rate  is  determined  by  a  complex  process  of  competition  between  the 
complete  set  of  transition  rates.  We  now  return  to  examine  the  photo¬ 
desorption  graphs  in  detail. 

The  high  temperature  behaviour  shown  in  Fig.  10  is  discussed 
first.  While  the  low  temperature  graphs  of  rate  versus  fluence  are 
qualitatively  similar  to  those  for  system  A,  the  desorption  rate  at 
T  >  150  K  and  low  fluence  is  limited  by  a  straight  diagonal  line.  This 
effect  is  a  characteristic  of  our  model,  and  is  easily  explained  using 
Fig.  12.  We  begin  with  the  adsorbate  molecule  in  the  ground  state 
j  =0,  v=0.  (For  temperatures  <150  K  the  initial  Boltzmann  energy  dis¬ 
tribution  dictates  that  practically  all  adsorbate  molecules  start  out 
in  the  ground  state.)  Since  the  lowest  surface  states  are  separated 
by  greater  than  one  Debye  energy,  our  one  phonon  transitions  are  cut 
off.  The  molecule  must  therefore  undergo  a  vibrational  (laser  induced) 
transition  to  the  j  =0,  v  = 1  state  before  proceeding  further  on  a  path 
to  desorption.  At  low  enough  fluences,  this  becomes  the  sole  rate 
determining  step,  governed  by  ~  10  F^.  The  diagonal  line  in 

Fig.  10  obeys  this  equation.  This  effect  is  manifested  in  Fig.  11  by 
the  rolling  off  of  the  rates  at  high  temperatures.  The  roll-off  points 
fall  on  a  straight  line  (dotted  line  in  Fig.  11)  that  may  be  charac¬ 
terized  by  an  activation  energy  of  6.9  This  is  very  nearly  the 

depth  of  the  j  =2  surface  state,  which  is  the  lowest  state  accessible 
via  one  phonon  processes.  The  dotted  line  appears  to  represent 
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thermal  desorption  from  the  j  =2  state!  With  this  observation,  we 
postulate  that  the  roll-off  occurs  when  the  vibrational  transiton  rate, 
L^,  is  exceeded  by  the  thermal  desorption  rate  from  the  j  =  2  state. 

At  high  temperatures,  the  main  path  to  desorption  apparently  involves 
three  general  steps,  as  shown  in  Fig.  13.  The  molecule  is  first  vib ra¬ 
tionally  excited,  then  couples,  via  phonon  emission,  to  an  excited 
surface  state,  and  finally  desorbs  via  a  thermal  cascade  of  phonon- 
assisted  steps.  The  roll-off  in  Fig.  11,  and  hence  the  corresponding 
diagonal  line  in  Fig.  10,  is  thus  merely  a  result  of  restricting  our 
theory  to  one-phonon  interactions,  while  the  lowest  surface  states  are 
spaced  by  greater  than  one  Debye  energy. 
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Fig.  13.  Path  to  desorption  at  high  temperatures 
for  system  B. 


At  lower  temperatures  now.  Fig.  11  shows  a  very  strong  enhancement 
of  desorption  rate  as  fluence  is  increased.  At  lower  fluence  levels, 
the  enhancement  is  less.  In  fact,  at  fluences  of  10  j/m  and  less, 
the  photodesorption  rate  at  T  <  75  K  drops  below  that  level  estimated 
several  paragraphs  back  by  ignoring  P  ,!  (The  solid  lines  in  Fig.  11 
fall  beneath  dashed  lines  for  F  <  10"^  J/m  .)  It  appears  that,  at 

AJ 

these  low  temperatures,  the  thermal  cascade  (see  left  side  of  Fig.  12) 

is  directed  strongly  downwards  because  spontaneous  emission  of  phonons 

dominates  in  the  rates  pTT, •  The  cascade  then  acts  as  a  sink  for 

13 

adsorbate  molecules,  resulting  in  decreased  desorption  rates  at  low  T. 
It  appears  that  the  temperature  dependence  of  the  desorption  rates  may 
be  attributed  mainly  to  the  characteristics  of  the  thermal  cascade. 

Finally,  at  the  lowest  temperature  -  T=  50 K  -  the  graphs  in  Fig. 11 
have  mostly  levelled  out,  to  become  independent  of  temperature.  This 
is  expected,  since  the  individual  transition  rates  gain  their  tempera¬ 
ture  dependence  via  the  phonon  density  of  states,  which  vanishes  at  low 
temperatures.  The  different  levels  for  different  fluences  at  T  =  50  K 
indicates  that  the  sink  mechanism  competes  less  efficiently  against  the 
desorption  processes  at  higher  fluences.  While  we  hope  that  the  simple 
interpretation  given  above  will  provide  some  qualitative  insight  into 
the  underlying  rate  mechanisms  involved,  it  is  emphasized  that  the 
actual  prediction  of  desorption  rates  necessitates  the  numerical  solu¬ 
tion  of  the  complete  master  equation.  This  completes  the  analysis  of 
system  B . 

The  photodesorption  rate  graphs  for  system  C  are  shown  in  Figs.  14 
and  15.  Since  the  bound  surface  states  of  this  system  are  very  widely 
spaced  (see  Fig.  7) ,  the  thermal  cascade  is  expected  to  play  a  minimal 
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Fig.  14.  Log  rate  versus  log  fluence  for  system  C 
at  various  temperatures. 
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Fig.  15. 


Log  rate  versus  inverse  temperature  for  system  C  at 
various  fluences.  (Units  of  fluence  are  j/m^.) 
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role.  As  a  result,  the  temperature  dependence  of  the  photodesorption 
rates  is  considerably  less  than  for  the  previous  two  systems. 

Comparing  the  three  systems  yields  important  qualitative  trends. 

As  y  ^  is  decreased  (from  system  A  to  B  to  C) ,  the  desorption  rate 
saturates  at  higher  levels  and  at  higher  fluence  values.  The  tempera¬ 
ture  dependence  decreases  as  y  ^  decreases.  The  number  of  bound  states 
also  decreases  with  y  \  eventually  preventing  thermal  desorption  from 
the  ground  state.  It  is  thus  observed  that,  by  varying  y  ^  among 
these  three  systems,  the  characteristic  fluence  and  temperature  depen¬ 
dence  of  desorption  is  dramatically  altered. 

Equally  importantly,  the  magnitude  of  the  desorption  rate  is 
changed  by  many  orders  of  magnitude  by  varying  y  \  In  Table  2  the 


Table  2. 

Desorption  rates 

(s  at  typical  temperatures 

and 

fluence s 

for  the  CH^F/NaCl 

systems  A,  B, 

and  C. 

System  (y  ^  in 

A) 

T(K) 

Fp ( J/m2) 

A( . 35) 

B ( . 25 ) 

C( .15) 

50 

.1 

1.4xl04 

4.5xl05 

9 

1.5x10 

50 

100 

1.5xl04 

2.9xl09 

7.6X1011 

150 

.1 

1.6xl07 

5.0xl08 

9 

1.5x10 

150 

100 

1.8xl07 

2.6xl010 

1.4xl012 

250 

.1 

9 

5.2x10 

9 

2.5x10 

9 

1.5x10 

* 
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desorption  rate  for  each  system  is  compared  at  several  values  of 

2 

fluence  and  temperature.  The  rate  at  F  =.l  j/m  and  T  =  50  K  de- 

X/ 

creases  by  5  orders  as  y  ^  is  increased  through  systems  C,  B,  and  A. 

-1  o 

If  we  increase  y  further  to  1  A  (then  =65.2;  r  =  440),  the 
fluence  dependence  disappears  and  photodesorption  is  no  longer  detec¬ 
table  above  the  thermal  rate,  even  at  T  =  50  K.  This  is  understood  to 

f  f  V 

result  from  the  fact  that  the  rates  P...,  P  and  Q  ..  become  very 

JJ  cjf  xcjf 

small,  exceedingly  quickly,  as  the  difference  between  the  states  - 

i  i  —1 

|j-jf|-  increases.  This  effect  is  amplified  as  oq  increases  with  y 

w 

As  a  result,  the  rates  P^+^  j  dominate,  inducing  step  by  step  transi¬ 
tions  in  the  thermal  desorption  cascade. 

The  correct  choice  for  y  is  thus  crucial  in  modelling  real  systems. 
Additional  spectroscopic  data,  such  as  bound  state  separation  energies, 
may  be  used  to  fix  y  more  accurately. 

We  note,  in  conclusion,  that  our  theory  depends  strongly  on  the 
shape  and  depth  of  the  surface  potential.  (While  the  adsorbate  molecu¬ 
lar  mass  (m)  and  the  Debye  temperature  (T^)  of  the  solid  are  also 
important  in  determining  and  r,  these  parameters  are  more  accurately 
known.  It  will  be  seen  -  through  examples  in  the  next  section  -  that 

the  size  of  the  fundamental  internal  molecular  vibrational  energy  (fioo  ) 

v 

in  comparison  with  the  well  depth  (V  )  is  also  important.  On  the  other 
hand,  the  other  characteristics  of  the  solid  and  the  adsorbate  enter 
the  rate  equations  in  a  much  more  innocuous  manner. 

§3.2  Comparison  with  Early  Experimental  Results;  Discussion 

We  will  now  compare  the  theoretical  data  of  the  last  section  with 
experimental  results.  Although  photodesorption  experiments  are  —  like 


. 


- 


60 


the  theories  -  in  the  earliest  stages,  the  kinetics  of  the  CH^F-on-Nad 
system  has  been  studied  by  Heidberg  et  al.  .  He  begins  with  an  adsor¬ 
bate  of  monolayer  to  multilayer  coverage  on  the  NaCl  surface,  at  64  K 
in  a  high  vacuum.  The  sample  is  exposed  to  pulsed  radiation  from  a 
high  power  TEA  C02  laser  which  has  been  tuned  to  the  vibrational  fre¬ 
quency  of  the  F-CH^  bond.  Desorbing  molecules  are  detected  with  a 
mass  spectrometer.  After  repeated  pulses,  the  adsorbate  coverage  is 
depleted.  At  a  pulse  fluence  of  the  order  of  F  ~  2-10  j/m  -  this 

X/ 

2 

is  equivalent  to  an  intensity  of  1  MW/ cm  -  he  measures  the  photo¬ 
desorption  rate  t  ^  z  10  s  with  the  fluence  dependence  given  by 


Our  model  does  not  immediately  predict  these  results.  A  compari¬ 
son  with  the  T=50  K  line  in  Figs.  8,  10,  and  14  shows  that  our  model 
predicts  too  large  rates  for  smaller  systems,  and  much  too  little 
fluence  dependence  for  the  larger  two  systems  which  are  in  saturation. 

The  maximum  fluence  dependence  observed  anywhere  in  our  systems  occurs 

2  -1  1.9 

in  system  B  at  T=50  K,  for  F^  <  1  J/m  ,  where  we  find  t^  “  F^'  . 

Noting  that  a  10%  increase  in  Vq  shifts  the  j  =0,  v=2  state  below 
the  continuum,  we  have  tested  the  effect  of  this  change  on  system  B. 
Increasing  to  25  kJ/mole  raises  oq  to  17.1.  We  now  take  vmax  =  3, 

and  then  find  that  the  desorption  rates  are  not  altered  appreciably! 

The  discrepancy  between  our  predictions  and  the  experimental 
results  does  not  necessarily  imply  an  inconsistency.  Our  theory  is 
designed  for  low  coverage  (that  is,  sub -mono layer)  systems,  in  which 
adsorbate  particles  do  not  interact  among  themselves.  The  experiments 
of  Heidberg  et  al.  involve  a  multilayer  adsorbate  which  will  likely  be 
influenced  by  interactions  between  adsorbate  molecules.  The  precise 


' 
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effect  of  such  adsorbate  molecule -molecule  interactions,  and  the  extent 
to  which  they  may  contribute  to  the  observed  discrepancy  are  not  known. 
Our  estimate  of  the  dipole  charge  separation  (  Q)  ,  taken  from  measure¬ 
ments  on  CH^F  in  the  gas  phase  as  previously  mentioned,  is  also  of 

questionable  validity.  Since  the  fluence  coordinate  scales  according 
2 

to  0  ,  an  error  in  Q  can  cause  a  significant  shift  along  the  fluence 
axis  in  our  graphs.  Such  a  shift  could,  for  example,  take  one  from  a 
region  of  fluence  saturation  to  one  of  considerable  fluence  dependence. 
This  remains  an  open  question  until  the  correct  value  of  Q  is  known. 

Another  feature  that  we  have  ignored  in  our  discussion  is  surface 
heating,  which  may  occur  in  two  ways.  First,  the  surface  may  absorb 
laser  radiation  directly;  this  effect  is  safely  neglected  in  an  essen¬ 
tially  lossless  dielectric,  such  as  NaCl  in  the  infrared.  Second,  the 
surface  absorbs  phonons  as  an  adsorbate  molecule  cascades  downwards 
through  the  surface  states.  The  sample  transition  path  shown  in  Fig.  16 
shows  how  it  may  be  possible  for  a  vib rationally  excited  molecule  to 
transfer  almost  all  its  energy  to  the  phonon  bath.  A  detailed  calcula¬ 
tion  of  the  surface  heating  via  this  type  of  mechanism  has  not  yet  been 
done . 

Phot ode sorption  experiments  on  somewhat  different  systems  have 
been  carried  out  by  Chuang  .  The  adsorbate  is  an  organic  ring  molecule, 
pyridine,  and  desorption  from  both  metal  (Ag)  and  dielectric  (KC1) 
surfaces  has  been  observed.  By  varying  the  laser  frequency,  they  find 
that  the  pho todesorp tion  yield  is  peaked  at  a  characteristic  vibrational 
mode  of  the  ring  molecule,  indicating  that  the  laser  couples  resonantly 
to  the  adsorbate  molecules.  On  the  metal  surface,  they  estimate  that 
direct  heating  of  the  solid  contributes  to  30%  of  the  desorption  at  a 
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Fig.  16.  A  laser-induced  heating  mechanism  for  system  B 


3  2 

fluence  F  =10  j/m  ,  and  to  only  10%  at  half  that  fluence;  the  contri- 

Xj 

bution  drops  off  rapidly  at  lower  fluences.  Chuang  concludes  that  the 

dominant  photodesorption  mechanism  here  must  involve  molecular  resonant 

absorption  of  the  laser  radiation.  We  do  not  attempt  to  apply  our  simple 

(diatomic  molecule)  theory  to  the  pyridine  system. 

Photo desorption  of  CO  from  a  Cu  crystal  by  non-resonant  surface 

4 

heating  has  been  recorded  by  Viswanathan  et  al.  Using  a  high  intensity 

3  2 

laser  with  a  pulse  fluence  F  s  10  J/m  ,  they  desorb  CO  molecules  at  a 

Xj 

-1  4  -1 

rate  t  ;  10  s  .  Our  theory  predicts  that  a  laser  with  pulse  fluence 
d 

2 

F  =1  J/m  tuned  to  the  CO  infrared  absorption  line  would  produce  about 
the  same  desorption  rate. (See  Section  3.3  for  details.)  This  prediction 


awaits  experimental  verification! 


’ 
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§3.3  Other  Theories 

It  is  worthwhile  to  compare  our  theory  presented  here  with  several 
others  in  the  literature. 

In  a  series  of  papers'^,  Lin  and  George  have  developed  a  theory  of 
photodesorption  based  on  a  reaction  mechanism  very  different  from 
ours.  They  suppose  that  the  laser  radiation  is  directly  absorbed  by 
the  surface  bond,  and  ignore  coupling  to  the  internal  molecular  vibra¬ 
tions.  (This  channel  is  shown  in  Fig.  lc  in  Section  1.)  This  energy 
absorption  then  induces  transitions  to  higher  energy  surface  states, 
with  the  assistance  of  the  phonons  of  the  solid.  The  surface  bond  is 
modelled  by  a  harmonic  potential,  plus  an  anharmonic  correction  term; 

desorption  is  defined  to  occur  when  a  molecule  reaches  a  preselected 

9 

energy  level  in  this  potential.  Our  criticism  of  this  model  is  two¬ 
fold: 

i)  Since  the  real  surface  potential  is,  in  general,  anharmonic, 
the  associated  energy  level  spacings  will  vary  considerably. 

The  laser  frequency  cannot  be  matched  to  all  these  spacings  at 
once!  The  laser  coupling  to  most  of  the  transitions  in  a  multi¬ 
state  potential  will  therefore  be  non- resonant ,  and  hence  quite 
weak.  (For  this  reason,  they  consider  also  a  two  laser  system, 

with  the  second  laser  tuned  to  excite  electronic  transitions  in 

5cs 

the  adsorbate  too!  ). 

ii)  The  validity  with  which  one  may  approximate  the  surface  bond 
with  their  very  simplified  model  potential  is  unclear,  and 
apparently  unnecessary,  as  the  following  authors  point  out. 

An  apparent  improvement  on  the  model  described  above  has  been  made 
by  Jedrzejek  et  al.6  They  consider  the  same  desorption  channel  as  Lin 


•  V  .. 
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and  George,  but  use  a  much  more  realistic  Morse  potential  to  model  the 
surface  bond.  However,  this  work  too  is  subject  to  critizism  (i) 
above.  As  a  result  of  weak  coupling  to  the  laser  radiation,  Jedrzejek 
finds  that  his  theory  needs  laser  intensities  four  orders  of  magnitude 
larger  than  realistic  values,  to  give  measurable  desorption  yields. 

In  contrast,  our  theory  predicts  desorption  rates  that  saturate  with 
respect  to  intensity  at  relatively  small  fluences.  The  predictions  of 
our  theory  may  be  compared  with  Jedrzejek’ s  in  Fig.  17,  where  data 
applied  to  the  CO-on-Cu  system  is  graphed.  Since  our  theory  predicts 
much  faster  rates  at  lower  fluences,  it  is  concluded  that  resonant  laser 
coupling  to  molecular  vibrations,  as  employed  in  our  theory,  initiates 
a  much  more  efficient  photodesorption  process. 

Photodesorption  rates  have  also  been  calculated  by  Lucas  and 
Ewing^,  employing  a  part  of  the  channel  used  in  our  theory.  They  only 
include  the  tunneling  term  -  essentially  -  between  excited  vibra¬ 
tional  states  and  continuum  surface  states;  see  Fig.  18.  For  the 
system  they  consider  the  vibrational  energy  spacing  is  much 

greater  than  the  surface  well  depth  (i.e.  heat  of  adsorption).  They 
assume  saturation  of  the  vibrational  states,  and  their  estimate  of  the 
desorption  rate  is,  in  their  approximations,  equivalent  to 


In  the  specific  regime  they  consider,  our  theory  reduces  to  approximate¬ 
ly  this  result  too.  This  simple  result  can  be  easily  understood  after 
noting  that  transitions  between  excited  vibrational  states  and  bound 
states  are  highly  unlikely  because  of  the  very  large  energy  difference. 
Consequently,  desorption  proceeds  mainly  by  the  two  parallel  processes: 


* 


log  rat 


log  F£  (J/m2) 


Fig.  17.  Log  rate  versus  log  fluence  for  the  CO/Cu  system. 
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Fig.  18.  Desorption  mechanism  for  the  H  /A1?0  system. 


i)  thermal  desorption,  and 

ii)  laser  assisted  via  (j,v)  =(0,0)  +(0,1)  +  (c,0). 

This  is  illustrated  in  Fig.  18.  At  the  low  temperatures  considered, 

(40  K)  thermal  desorption  is  negligible,  so  process  (ii)  dominates.  For 

high  enough  laser  fluence,  the  vibrational  transitions  are  very  fast, 

and  hence  the  final  jump  to  the  continuum  is  the  single  rate-de te mining 

step.  In  our  theory,  the  rate  of  this  transition  is  just  +  P^q. 

Normally,  <  0^;  and  we  thus  get  almost  the  same  result  as  (97). 

cU i) ii) * *  vcU 

9 

Finally,  we  note  that  the  authors  Kreuzer  and  Lowy  have  previous¬ 
ly  considered,  instead  of  the  0^  term,  the  effect  of  the  tem. 

This  work,  which  involved  some  calculational  simplifications,  is  in  fact 
the  beginning  out  of  which  the  theory  presented  here  has  evolved. 
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§4.  CONCLUSIONS 


A  fully  quantum  mechanical  theory  of  desorption  via  laser- 
stimulated,  resonant  excitation  of  adsorbate  molecular  vibration  has 
been  presented.  The  theory  is  based  on  a  simple  microscopic  model  in 
which  desorption  entails  a  cascade  of  surface  and  vibrational  state 
transitions.  All  physical  parameters  introduced  in  our  model  may  in 
principle  be  fixed  by  experiments;  desorption  times  are  predicted  over 
a  full  range  of  temperature  and  fluence.  A  comparison  with  other 
theoretical  calculations  and  experimental  results  indicates  that  the 
resonant  laser  coupling  process  considered  here  may  be  a  useful  way  to 
selectively  and  efficiently  stimulate  desorption.  In  the  cases  cited 
in  Section  3,  the  use  of  a  laser  to  excite  internal  molecular  states 
appears  to  be  much  more  practical  than  exciting  the  surface  states 
di  re  ct  ly . 

An  actual  verification  of  the  present  theories  awaits  additional 
experiments;  obtaining  accurate  desorption  measurements  is  clearly  not 
an  easy  task.  Our  theory  relies  on  experimentally  determined  parame¬ 
ters,  such  as  T  ,  V  ,  id  ,  F  .  Q,  and  Y;  the  last  two  of  these  are  not 

D  o  v  Z 

presently  available  for  most  systems.  To  check  our  theory,  future 
photodesorption  experiments  should  be  performed  at  measured,  submono¬ 
layer  coverage.  The  potential  range  parameter  (y)  may  be  determined 
either  spectroscopically  or  by  thermal  desorption  experiments  done  at 
the  same  coverage.  (The  latter  method  provides  an  evaluation  of 
too.)  The  molecular  dipole  separation  (0)  may  in  principle  be  found  by 
measurement  of  the  infrared  absorption  intensities  of  the  adsorbate. 

An  unequivocal,  experimental  evaluation  of  our  theory  for  resonant 
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photodesorption  could  then  be  obtained. 

Several  necessary  improvements  to  our  theory  are  immediately 
apparent.  Most  notable  is  our  neglect  of  the  surface  heating  effect 
due  to  phonon  emission  from  inelastic  transition  processes.  The 
seriousness  of  this  omission  should  be  determined  by  future  calcula¬ 
tions.  Another  weakness  is  our  use  of  a  one  dimensional  surface 
potential  and  a  fixed  molecular  orientation.  A  complete  three  dimen¬ 
sional  calculation  (permitting  molecular  rotation)  seems  appropriate 
here.  It  must  also  be  remembered  that  our  theory  is  designed  for 
localized,  physisorbed  (weakly  coupled)  systems  at  negligible  surface 
coverage.  In  practice,  the  coverage  is  usually  greater,  allowing  for 
significant  interaction  between  adsorbate  molecules.  Resonant  "V-V" 

coupling  between  molecular  vibrational  states  in  the  adsorbate,  as 

7c 

suggested  by  Heidberg  for  the  CH^F/NaCl  system,  is  an  important  fea¬ 
ture  that  has  been  ignored  in  the  work  presented  here.  Furthermore,  we 
have  treated  the  solid  as  a  continuum  throughout,  and  coupled  the 
macroscopic  electric  field  of  the  laser  to  the  molecular  vibrations  of 

the  adsorbate  via  its  infrared  absorption  spectrum.  The  currently 

22 

developing  topic  of  "Surface  Enhanced  Raman  Spectra",  although  ignored 
here,  may  provide  additional  features  to  be  utilized  in  future  photo¬ 
desorption  processes  on  metal  surfaces.  There  are  clearly  many  areas 
open  to  study,  both  theoretical  and  experimental. 

It  summarizes  to  say  that  the  utility  of  photodesorption  -  via 
resonant  laser-molecular  coupling  -  as  a  tool  for  selective  desorption 
awaits  the  developments  of  many  ongoing  investigations. 
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APPENDIX  A 
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Orthogonality  of  (U_*  } 

Equation  (14)  is  now  derived  for  the  modes  (x)  defined  in 

ka 

equations  (13a, b) .  Using  (13a), 


dx  n  (x)U->  (x)*U-*.  (x)  = — ~ 


k3 


k '  3 '  4tt 
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dxn  (x)fA.  (x)*f4.  (x) 

KfCf 


dx  n  (x)l^.  (x)-f+  (x) 

KC  KC 


(Al) 


The  remaining  integral  in  (Al)  must  be  evaluated  explicitly  for  the 


cases  3  =  s,p.  Using  equation  (13b)  in  the  s-polarization  case  gives 
00  0 
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where  a'  and  b1  are  evaluated  at  £*.  To  arrive  at  the  second  line 
above,  we  note  that 


i(4’-C)x 
dx  e  =i 


r  +  ir6(c,-0 


0 


where  (p  denotes  the  principle  part  distribution.  The  terms  involv¬ 
ing  (r/C-41  all  cancel  in  equation  (A2)  .  Since  we  require  4,4’  <0,  the 

s 

term  in  (A2)  proportional  to  <5(c’+c)  does  not  contribute.  With  N 
defined  as 

1  +  n  a 

NS  =  r  -  (A3) 

/2tt(1  +2n  a  +  n2a2) 

r  r 


71 


equation  (A2)  reduces  to 


dx  n2(x)  fl  (x).fj.(x)  =  «(?’-?) 


(A4) 


Similarly,  substitution  of  (13b)  into  (Al)  for  the  p-polarization  case 
and  performing  analogous  manipulations  yields 


dx  n2(x)  f$  (x)-fj  (x)  =  6(c'-C) 
KC  K£ 


(A5) 


where 
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/2tt(K^  +  £2)  +  2n^a  +  a^) 


(A6) 


Finally,  substituting  (A4)  and  (A6)  into  (Al)  gives  the  orthogonality 
relation 


dx  n2(x)Uv  (x)U*  (x) =  6  f S(K-Kr )<5 (<;-£' ) = 5  , 6 (k-k' )  .  (A7) 

k3  k  3'  pp  pp 

Equation  (15)  is  derived  in  a  similar  manner.  Using  (13a),  we  have 
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Evaluating  the  integral  over  x  using  the  same  method  as  before  leads  to 


dx  n2(x)U^.  (x)*U^_  (x)  =  -  A_^  6  (K+Kf )  6  (  £-  C  f ) 
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where 


Aits 


.  2  2 

1  -  n  a 
r 


2  2 

1  +  2n  a  +  n  a 
r  r 


and  A->  = 


2  2 
n  -  a 
r 


k-P  n2  +  2n  a  +  a2 


(A10) 


. 


xh 


APPENDIX  B 


The  Morse  Wave  Functions 

The  eigen-solutions  to  the  Schrodinger  equation  for  the  Morse 


poten tial, 
*2  , 


-2y(x-x  ) 


2m  ^$(x)  +  Ve 


-y(x-x  ) 

2e  °  )$(x)  =  E$(x) 


(Bl) 


are  quoted  directly  from  reference  1.  The  energy  eigenvalues  take  on 
discrete,  bound  state  levels  given  by 


2  2 

E  =  JLX-s2 
J  2m  j 

,  .  1  .  n  i  12 
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(B3) 


where  q  is  the  wave  vector  of  the  particle. 


The  bound  state  wave  functions  are 


$  .  (x)  =  /y  <j> .  (c) 
J  J 


(B4) 


where  C  =  y(x-x  )  and 
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(B5) 


in  terms  of  the  Laguerre  polynomial  L_.  ^(u). 

The  continuum  state  wave fun c ti ons ,  normalized  in  a  box  of  length 
L,  are  given  by 


*  (x)  =  L  2  4^(5) 


where  n  =  q/y  and 
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with  the  confluent  hyper  geometric  function  V(a,b,z)  defined  as 

1  iTT  fFC^-o  +in,l+2in,z) 

¥(:«-- cr  +in,l+2in,z)  =  ■.  — r  1  w'-T.o-  \r/i - : — r  - 
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-2in  F(^"a0-ih»1-2in,z) 
z  r(i-2in)r(^-aQ+in)  j 

F(z,b,z)  is  a  Rummer’s  function,  as  defined  in  reference  23c 
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APPENDIX  C 


Evaluation  of  the  Integral  I 

^  J  3 

The  integral  I  ,  may  be  rewritten  by  substituting  the  wave 

J 

functions  (B5)  into  (76b)  to  get 
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(We  have  substituted  u  =  2a  e  ^.) 

o 

r00 

The  method  of  evaluation  of  )  ,  described  below  follows  that 

used  by  R.  Teshima  to  do  similar  integral  in  reference  1.  We  first 
separate  (C2)  into  two  parts 
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Since  I  is  symmetric  with  respect  to  j,j’,  we  let  j  >j*.  We  then 
°  2aQr2 j-1 

expand  L^.  (u)  using 
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and  then  integrate  Iq  by  parts  j  times  to  get 
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Similar  operations  performed  on  1^  lead  to 
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Now,  substitution  of  (C6)  and  (C7)  into  equation  (C4)  gives 
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We  finally  get  I  ,  by  substitution  into  equation  (Cl). 

3 

The  corresponding  integral  including  bound  and  continuum  state 
wave  functions  (I  is  almost  identical  to  one  appearing  in 

reference  1,  so  we  just  quote  the  final  result  here: 
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where  q  =  q/y. 
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